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This colloquium gives an overview of recent theoretical and experimental progress in the area 
of nonequilibrium dynamics of isolated quantum systems. We particularly focus on quantum 
quenches: the temporal evolution following a sudden or slow change of the coupling constants of 
the system Hamiltonian. We discuss several aspects of the slow dynamics in driven systems and 
emphasize the universality of such dynamics in gapless systems with specific focus on dynamics 
near continuous quantum phase transitions. We also review recent progress on understanding 
thermalization in closed systems through the eigenstate thermalization hypothesis and discuss 
relaxation in integrable systems. Finally we overview key experiments probing quantum dynamics 
in cold atom systems and put them in the context of our current theoretical understanding. 
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I. INTRODUCTION 

In the past two decades the outlook of condensed mat- 
ter physics has been deeply and unexpectedly revolu- 
tionized by a few experimental breakthroughs in atomic 
physics, quantum optics and nanoscience. In synthesis, 
crucial advances in these fields have made it possible 
to realize artificial systems (e.g. optical lattices, quan- 
tum dots, nanowires) that are described to a very high 
degree of accuracy by models (e.g. Hubbard, Kondo, 
and Luttinger models) whose physics has been a subject 
of intense investigation in various contexts ranging from 
high temperature superconductivity to low temperature 
transport in metals. It is fair to say that this experi- 
mental progress has changed the way theory and exper- 



iment look at each other. In the past, effective models 
were largely devised to explain the low energy physics 
of highly complex systems. The situation has now been 
reversed so that one can experimentally realize and simu- 
late the physics of such models. On one hand, the design 
and realization of interacting many-body systems could 
in principle be used to perform practical tasks, su ch as 
quantum information processing (Farhi et al. 2001). On 



the other hand, direct simulations of simple models could 
help resolving important problems in condensed matter 
physics. But most importantly, the availability of ex- 
perimental controllable systems whose properties can be 
accurately described by simple models provides unprece- 
dented opportunity to explore several new frontiers of 
condensed matter physics including the nonequilibrium 
dynamics in closed interacting quantum systems. 

Equilibrium systems can often be understood using 
a combination of a mean field theory, renormalization 
group, and universality. This allows us to understand 
low temperature experimental data obtained in complex 
systems, such as interacting electrons in solids, in terms 
of simple effective models containing a few relevant pa- 
rameters. Away from equilibrium the situation is much 
less clear. While some progress was made in the past for 
classical systems (Schmittmann and Zia 1995 1, there are 



no rigorously justified generalizations of any of these ap- 
proaches to generic quantum nonequilibrium systems. It 
is thus not obvious that the theoretical study of the dy- 
namics of simplified models would accurately describe ex- 
periments of more complex systems. In addition there are 
fewer available tools for analyzing dynamics of even sim- 
ple interacting models. In this respect cold atomic gases 
and nanostructures make possible what would be arduous 
otherwise: a fruitful comparison between nonequilibrium 
theories based on simple models and carefully designed 
experiments with tunable system parameters. 

Finding systematic ways to understand the nonequilib- 
rium physics of interacting systems is not only of funda- 
mental importance, but could also be crucial for future 
technologies. A quantum computer, for example, will 



definitely require the capability of performing real time 
manipulations of interacting quantum systems. Though 
large scale quantum computers are yet to be on the hori- 
zon, it is evident that quantum coherent dynamics will be 
one of the focus points of various experimental systems 
and of future technologies. 

Nonequilibrium dynamics is a potentially a vast field: 
there are many ways to take a system out of equilibrium, 
such as applying a driving field or pumping energy or 
particles in the system through external reservoirs as in 
transport problems. It is thus of utter importance to 
focus on simple, yet fundamental protocols. In this Col- 
loquium we will concentrate on the simplest paradigm: 
the study of the nonequilibrium dynamics of closed inter- 
acting quantum systems following a change in one of the 
system parameters (quantum quench). Such a change, 
which could be either fast or slow, is particularly in- 
teresting when it takes the system through a quantum 
phase transition involving macroscopic changes in the 
state of the many-body system at the initial and the final 
point. Seminal work in this direction includes ground- 
breaking experiments ( [Greiner et al. 2002a|b ) showing 
both the feasibility of observing a quantum phase transi- 
tion in cold atoms and the possibility of observing quan- 
tum coherent dynamics. Following this work, a number 
of different experiments explored the dynamics of cold 
atom systems driven across BCS-BEC crossover (Re- 



gal 2006 ) , polar and ferromag netic phases of spinor con- 
densates dSadler et a/.[|2006[), insulating and superfluid 



phases of ultracold bosons ([Tuchman et al. 2006 ) and 



many others (see Ref. ( Bloch et al. 2008 ) for a review) 



These experiments stimulated an active theoretical re- 
search in the relatively unexplored area of quantum dy- 
namics in closed interacting systems. An interesting 
characteristic common to these systems is that despite 
of the absence of energy exchange with an environment 
and of the consequent global relaxation, it neverthe- 
less frequently possible to look at the long time dynam- 
ics and characterize it in terms of an asymptotic state 
attained by physical (measurable) observables (Cramer 
'erar|[2008]|Flesch et aZ.j,2008,|Gogolin et aL[|20il[|Lin 



den et aL[ [2009[ |Reimann[ |2008[ |Rigol et al.\\'m^ . In 
connection to this, it is possible to categorize recent re- 
search on the subject of this Colloquium in two main 
questions: 

• What is universal in the dynamics of a system fol- 
lowing a quantum quench ? 

• What are the characteristics of the asymptotic, 
steady state reached after a quench ? When is it 
thermal ? 

In this Colloquium we will discuss both of these ques- 
tions extensively. We shall outline our current level of 
understanding of these issues and chart out the outstand- 
ing open questions in the field. In Sec. |ll] we will focus 
on the first question and describe, from various points of 
view, the universal aspects of nearly adiabatic dynamics 



near quantum critical points as well as in generic gapped 
and gapless systems. We will argue that the proximity 
to the adiabatic limit allows us to make specific univer- 
sal predictions of scaling of various quantities such as the 
defect density and heat with the quench rate. 

In Sec. |III[ we will discuss recent progress in under- 
standing thermalization of a quantum system following a 
quench. In classical systems active interest in this topic 
was stimulated by the celebrated work of Fermi Pasta 
and Ulam on the dynamics of a one-dimensional (ID) 
anharmonic chain (Campbell et al. 2005 Fermi et al.\ 
May 1955 Porter et al. 2009) which demonstrated the 
absence of such thermalization. It was realized much 
later that the nonlinearity of the interaction is not suf- 
ficient for thermalization which occurs, in this system, 
only if the initial amplitude of interaction exceeds a cer- 
tain threshold (Izrailev and Chirikov |1966 ). Below this 
threshold, the solution splits into solitons and retains 
its quasi-periodic nature ( [Zabusky and Kruskal[ 1965) 
which is a consequence of the Kolmogorov-Arnold-Moser 
(KAM) theorem (Tabor 19891. In quantum systems. 



the question of sufficient criteria for thermalization has 
remained largely unaddressed so far. Some experimental 
progress in this direction has been made by a recent ex- 
periment from Kinoshita et. al. ( Kinoshita et al. 2006 ) 



on non-thermalizing dynamics of ID bosons with short 
range interactions. This experiment constitutes the first 
clear demonstration of the fact that a nearly integrable 
quantum interacting many-particle system does not ther- 
malize for a very long time. Currently, the question of 
extension of KAM theorem to quantum systems is a sub- 
ject of active theore tical debate (see e.g. Ref. (Olshanii 
and Yurovsk^|2009[ )). 

Finally, we note that many important topics concern- 
ing the physics of closed interacting systems did not 
find space in this Colloquium. Most important among 
these are the tools that are being developed to de- 
scribe theoretically the physics of interacting systems 
out of equilibrium. Among such methods we men- 
tion density-matrix rcnormalization group (DMRG) and 
time-evolving block decimation (TEBD) for analyzing 



equilibrium and nonequilibrium ID systems (jSchollwock 
20051 ISchollwock and Whitel [2006HVidaI| 



White 



et al. 



1992) and higher dimensional ones (Verstraete 
2008), the Keldysh technique (Kamenev and 



2003 2004 



Levchenkol 2009 ) which is particularly helpful for deriv- 



ing quantum kinetic equations, and closely related func- 



tional integral methods (Gasenzer 2009 Plimak et al. 



2001 Rey et al. 2005 ) . Cold atom experiments also 



prompted rapid developments in phase space methods, 
where quantum dynamics is represented as an evolu- 



tion in the classical phase space (Blakie et al. [2008 
2010). These methods were originally de- 
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veloped and applied to various problems in single-particle 



dynamics (Hillery et al. 1984 Zurek 2003) and inde- 



pendently in quantum optics in the context of coherent 



states (Gardiner and Zoller 2004 Walls and Milburn 



1994). There are other reviews available in literature 



(see the references above) which specifically target these 
areas. 



II. NEARLY ADIABATIC DYNAMICS IN QUANTUM 
SYSTEMS 

A. Universality in a nutshell 

Universality (or insensitivity to microscopic details) is 
one of the crucial concepts of modern condensed matter 
physics. It naturally emerged from one of the milestones 
of modern physics: the renormalization group. In con- 
densed matter physics universality is a very powerful tool 
for the understanding of co ntinuous (second order) phase 
transitions, both classical ( Chaikin and Lubensky' '1995', 



Landau and Lifshitz. 1980 



Sondhi et al. 



1997 Vojta 



) and quantum ( Sachdev,|1999} 
2003 1 . As a consequence of the 



divergence of the correlation length, a system undergo- 
ing such a continuous phase transition is typically scale 
invariant in the vicinity of the critical point and can be 
characterized by relatively simple massless field theories, 
which permit a classification of perturbations driving the 
system away from the critical point. Consequently, uni- 
versality manifests itself in the scaling behavior of var- 
ious quantities such as the order parameter, (free) en- 
ergy, susceptibilities and correlations functions near the 
critical point. In this review we will focus mostly on 
quantum phase transitions occurring at zero tempera- 
ture upon the variation of a control parameter A through 
a critical point Ac. A standard example of universality 
is the fact that the exponent i/ characterizing the diver- 
gence of the correlation length ^^1/|A — Ad'' near the 
quantum critical point (QCP) is insensitive to the mi- 
croscopic details of the system and depends only on the 
universality class of the transition, determined by the di- 
mensionality, overall symmetries and range of the inter- 
actions. For classical (thermal) phase transitions similar 
universality manifests in the divergence of the relaxation 
time Trei ^ 1/|A — Acl^'^, where z is the dynamical crit- 
ical exponent. For quantum phase transitions the expo- 
nent z can be associated with a vanishing energy scale 
A ~ |A — Acl^", which can be either a gap or a crossover 
scale where the spectrum changes qualitatively. By the 
uncertainty principle this energy scale corresponds to 
a divergent time scale, which typically describes the 
crossover in the scaling behavior of unequal time cor- 
relation functions. Phase transitions can be also char- 
acterized by singular susceptibilities, which are in turn 
connected through the fluctuation-dissipation theorem to 
the correlation functions of conjugate variables (e.g., the 
magnetic susceptibility is related to the correlation func- 
tion of the magnetization) . At critical points these corre- 
lation functions have often power law scaling behavior at 
long distances, e.g. {m{x)m{x')) « l/\x — x'\'^". The ex- 
ponent a sets the scaling dimension of the corresponding 
operator m(x): d[ia[m{x)] = a. Because similar correla- 
tion functions can enter different susceptibilities not all 



the scaling exponents are independent but must satisfy 
scaling relations (Chaikin and Lubensky 1995 Vojta 
2003|. 



As mentioned in the introduction, the idea of univer- 
sality makes it possible to interpret experimental data ob- 
tained in real systems in terms of effective models with 
a few parameters. Universality can be ultimately un- 
derstood using the renormalization group, which shows 
that as a system is coarse grained to lower energies and 
longer length scales, more and more parameters of its 
original, ab initio description become unimportant (ir- 
relevant), while the remaining few (relevant) parameters 
define an effective low energy model. A standard example 
of universality in this context is the scaling relation be- 
tween energy and momentum of elementary excitations, 
e oc k^ , controlled by the dynamical exponent z which 
depends on the symmetries of the system. In particular, 
z = 1 in most phases characterized by a continuous bro- 
ken symmetry (crystals, superfluids, anti-ferromagnets) , 
z = 2 in ferromagnets, where there is an additional con- 
servation law of the order parameter. 

Universality is well established and understood in equi- 
librium. It is, however, crucial for many experimentally 
relevant situations to understand the extent to which this 
concept can be extended to out of equilibrium physics. 
Can irrelevant interactions turn out to be important away 
from equilibrium? Since there are many ways to take a 
system out of equilibrium, for which specific protocols 
will universality emerge and which details of the proto- 
col are potentially important? Below we will focus on re- 
cent studies addressing these important issues in closed 
interacting quantum systems, and in particular on the 
dynamics of a system whose parameters are dynamically 
tuned either through a quantum critical point, or in gen- 
eral within a gapless/gapped phase. 



B. Universal dynamics near quantum critical points 

Let us start by considering the simplest nonequilibrium 



protocol (Dziarmaga 2005 Polkovnikov 2005 



et al. 2005): the system is prepared in its ground state 



Zurek 



and is then driven through a QCP by changing an exter- 
nal parameter A in time. As long as the rate change of the 
gap in the spectrum A ^ caused by changing A is smaller 
than the square of the gap one can expect the system to 
approximately follow the ground state adiabatically (we 
will revisit this statement in the next section). However, 
the vanishing of the gap at A = Ac implies that the sys- 
tem will always violate adiabaticity close to the quantum 
critical point, no matter how slowly the parameters are 
changed. It is then natural to ask how many excitations 



^ We use the word "gap" for brevity. However, the system can be 
gapless on one or both sides of the transition (e.g. superfluid- 
insulator transition). Then A would denote a crossover energy 
scale vanishing at the QCP. 



will be generated while passing though the critical point 
and how their density as well as generated entropy and 
energy will depend on the rate of change of A. 

A similar question known under the name of the 
Kibble-Zurek mechanism (KZ) 



Time after quench 



(Kibble 2007 1976 



Zurek 1985 1996 1 has been addressed in the past decades 
for classical phase transitions. In that case the excita- 
tion density of defects is described by a simple scaling 

Suppose that the tuning 



1985 1996) 



argument (Zurek 
parameter, for example external temperature T, slowly 
decreases in time across the critical value T^- T = Tc — vt. 
The system will respond adiabatically (quasi-statically if 
the system is not thermally insulated) up to some close 
vicinity of the critical point, where adiabaticity will be 
violated as a result of the divergence of the relaxation 
time {Trei ~ ^/\T — T^Y'^) and the dynamics will become 
diabatic (sudden). The adiabatic response is once again 
resumed after the system moves out of the vicinity of 
the critical point. Zurek suggested a very simple crite- 
rion for separating such adiabatic and diabatic (impulse) 
regimes: the time to reach the critical point t = \T—Tc\/v 
should be equal to the relaxation time. This immedi- 
ately introduces the time and length scales characterizing 
the adiabatic to diabatic crossover: t* ~ 1/|-l'|^'^/^^'^+^', 
S,* ~ l/|-u|'^/(^''+^^. The violation of adiabaticity implies 
that order can not form on distances larger than ^* lead- 
ing to the formation of a domain structure with a char- 
acteristic distance ^* between the domain boundaries. 
In two and three dimensional systems, when the order 
parameter is characterized by a continuous broken sym- 
metry, the points where several domains meet correspond 
to vortices or vortex rings. These are robust topological 
excitations with a very long life time (see Fig. [I]). Since 
^* determines the average distance between the defects 
their density is given by a simple universal expression 



*\d 



\di'/(zi^+l) 



(1) 



The universality of the KZ prediction above is manifest in 
the appearance of the universal critical exponents z and i^ 
in the scaling law. This scaling was confirmed in experi- 
ments in liquid crystals ( Ducci et al. |1999 ) . Experiments 
in other systems (superconductors (Maniv et al. 2003), 



arrays of Josephson junctions (Monaco et al. 2006)) ob 



served the production of topological defects with a power 
law scaling on the quench rate but gave a different ex- 
ponent. The KZ scaling was also confirmed theoreti- 
cally using stochastic dynamics (Ginzburg-Landau dy- 
namics with Langevin noise or Glauber dynamics) where 



temperature changes in time (Krapivsky, 2010; Laguna 
[and Zurekl |1997[ [Rivers} |200l| [Yates and Zurek^ ,1998) , 




FIG. 1 Defect generation after a quench in a spinor 
condensate - These images show the transverse magnetiza- 
tion density of spinor condensates for variable evolution times 
after a quench to a ferromagnetic state, revealing a spatially 
inhomogeneous formation of ferromagnetic domains. The ori- 
entation <j) and amplitude A are depicted by the color and 
brightness according to the color wheel shown. Inset: An 
instance of a spin vortex spontaneously created during the 
quench. For reference, the length scale corresponding to the 
characteristic healing length ^ is also shown. (Adapted from 
Ref. (iSadler et aZ.||2006|). See Sec. 



IV for more details. 



Finding the scaling of these quantities remains an open 
question. 

The arguments above were recently generalized to the 



though there are also works suggesting various modi- crossing of quantum phase transitions (Dziarmaga 2005 



fications (Biroli et al. 2010 Hindrnarsh and Rajantie Polkovnikov 2005 Zurek et al. 2005 1 (for recent reviews 



2000). One can also interpret Eq.([lJ) as a measure of 
non- adiabaticity near the critical point. It is naturally 
expected that other measures like non-adiabatic energy 
production and entropy generation will display similar 
universality. These measures might be preferable over 
ricx in situations where it is difficult to identify defects. 



on this subject see Refs. (Dziarmaga 2010 Gritsev and 



Polkovnikov 2010). As discussed before in the quantum 



case the parameter to be varied is not temperature T 
but rather the coupling A tuning the system through the 
quantum critical point. In order to obtain the scaling 
for the number of excitations produced in the quantum 



case let us first recall the Landau-Zener analysis of the 
crossover between adiabatic and nonadiabatic dynamics 
in a simple driven two-level system: 



'H-lz = 9{t)(^z + Acr:r-, 



where g{t) — vt. If the system was initially prepared in 

the ground state at f — > — oo, the proba bility of t ransi 

tion to the excited state at i — >■ +00 is (Landau 1932 



Majorana[ [l932[ |Stiickelberg[ |1932( |Zener| |1932f 

Pex = exp[-7r7], (3) 

where we introduced the Landau-Zener parameter 7 = 
A'^/v. Notice that the limit 7^1 corresponds to the 
adiabatic limit with an exponentially suppressed tran- 
sition probability while 7^1 corresponds to the dia- 
batic limit where the transition happens with probabil- 
ity close to unity. Hence when the rate of change of 
the energy splitting between two levels becomes larger or 
comparable to the energy splitting squared one observes 
a crossover from adiabatic to diabatic dynamics. An al- 
ternative qualitative explanation of this result has bee 
recently formulated (Damski and Zurek 2006). 



The Landau-Zener argument can be straightforwardly 
extended to the crossing of a QCP. The characteristic 
energy scale which changes in time is now the gap A. 
As we discussed earlier this gap universally depends on 
the tuning parameter A: A(A) ^ |A — Xc\^'^ ~ |wi|^'^, 
where we assumed that the dependence X{t) can be lin- 
earized near the QCP: A(i) ?« Ac -I- vt. Comparing the 
rate of change of the gap with its square, i.e. solving 
the equation dA/dt sa A^, we find the energy scale at 
which adiabaticity breaks down is A* ^ |'li|^'^/(^''+i). At 
this point the system is characterized by the length scale 
C* ~ \v\~'^^^'"^^^\ which can be interpreted as the typ- 
ical length scale of fluctuations of the order parameter. 
Beyond this point the adiabatic approximation breaks 
down and fluctuations at longer scales cannot adiabati- 
cally follow the ground state. This results in the creation 
of defects with typical distance ^ between them and den- 
sity nex ~ Id'' ~ \v\d''/i^'^+^). This scaling is identical 
to the classical one, Eq. Q with A — » T, and was pro- 
posed indep endently in Refs. ( [Polkovnikov 2005 Zurek 



et al. 



2005 ) . There is a simple quasi-particle interpreta- 



tion for this scaling: assuming that the excitations in the 
system are characterized by isolated quasi-particles then 

their density can be found from riox ~ Jr. dep{e), where 
p(e) is the single-particle density of states near the QCP. 
In uniform d-dimensional systems p{e) ~ ^d/z-i^ which 
again reproduces Eq. pj). 

The scaling in Eq.(fll was verified in a series of ex- 
act solutions of the dynamics across the QCP in in- 
tegrable systems whose dynamics can be mapped into 
a series of Landau-Zener transitions of a few quasi- 
particle modes. In particular, it has been verified for 
various spin models in one and two dimensions which 
can be mapped to noninteracting fermions (Mukherjee 
et al. 2007 ) , for models where low energy excitations 



near phase transitions can be described by bosonic Gold- 



stone modes (Dziarmaga et al. 2008 Lamacraft 2007 



Polkovnikov 2005 1 , the sine-Gordon model, where el- 
ementary excitations are solitons and breathers with 



(2) fractional statistics (iDe Grandi et all 120081 |2010b|). 



graphene (Dora et al. 



2010 Dora and Moessner 20101. 



This scaling was also extended to disordered systems, 
like a disordered Ising spin chain, where it was found 



that riex ^ 1/log (u) (Caneva et al. 



2007 



Dziarmaga 



2006[ ), as expected from Eq. (|1| due to the divergen ce of 
the exponent z near the critical point (Fisher 



1995). 



The scaling in Eq. (fTl) can be generalized to the case of 
nonlinear dependence of the tuning parameter on time, 
A(i) ^ Ac ± ^^1^, whe re considerations similar to those 
leading to Eq. ([l|) give ( B arankov and Pol kovnikov [2008 



De Grandi et a?.[[20T0a[|Sen et a/.[|2008p : 



\di'/{zi/r-\-l) 



(4) 



In all cases v in Eq. ^ plays the role of the adiabatic 
parameter: the limit f — >■ corresponds to the adiabatic 
limit (this interpretation is valid even for instantaneous 
quenches r = 0, where v plays the role of the quench 
amplitude). This suggests the dynamics can be sys- 
tematically analyzed using adiabatic perturbation theory 



(De Grandi and Polkovnikov 2010 Polkovnikov 2005 



Rigolin et al. 2008), i.e. expanding the transition am- 



plitudes to the instantaneous eigenstates of the system 
in powers of v. Using such analysis in Refs. (De Grandi 
et al. 2010a|b ) it was shown that the scaling (|4| can be 
derived from the scaling of the adiabatic fidelity, defined 
as the overlap of the time-dependent wave function with 
the instantaneous ground state: 



F(i) = |(V'(t)IV'gsW)l- 



(5) 



In particular, for A(i) = Ac -I- vV /r\0{t)^ where 0{t) is a 
step function, 



P^V) = 1 - F{tY « L''\v\'x2r+2{\c), (6) 



where 



X2r+2(A) 



|(n|^|0)|^ 



U^ ^-^ (En 



^0) 



2r+2 ■ 



(7) 



is the adiabatic (fidelity) susceptibility of the order 2r-|-2 



(X2 is the conventional fidelity susceptibility ( Gu and Lin 
2009)). Here En are the eigenenergies and V is the op- 
erator coupled to the parameter X: V = d\H{X)\x=\^. 
If the perturbation is local and spatially uniform , i.e. 
V = J d'^xu{x), then the scaling dimension of the adia- 
batic fidelity susceptibility is obtained from a straight- 
forward generalization of the result of Refs. (Gu and 



2009[|Venuti and Zanardi| |2007[ ) , i.e. dim[x 



<.2r+2i 



2z(r -I- 1) — d, where A„ = dim[M] is the scaling 



Lin 

dimension of u(a;). 

Let us now discuss from this general perspective the 
arguments leading to the generalized scaling relation 



Eq. Q . If the scaling dimension of the susceptibihty is 
negative, this imphes that X2r+2 diverges at the critical 
point. In this case from Eq. m we find that asymptoti- 
cally at w — ^ 



PcAy) 



|2 r2d+2z(r+l)-2A„ 



(8) 



From Eq. ([8]) we see that the probability of excit- 
ing the system becomes of order one when L ^ 
l/\y\i/id+z{r+i)-A^)^ This length scale can be inter- 
preted as the typical distance between elementary exci- 
tations (defects) and thus we find that instead of Eq. (Is]) 
we get 



|d/(d+2(r+l)-A„) 



(9) 



This expression reduces to Eq. (HI) if u(x) is a relevant 
operator driving the system to the new phase. Indeed, 
in this case X J d'^xu{x) should have the same scaling di- 
mension as the gap, i.e. z, which immediately implies 
that the scaling dimension of u{x) is Am = d + z — \/v 
and that dim[x2r-i-2l = d— 2z r — 2/v (|De Grandi et al. 



2010b Schwandt et al. 2009). Notice finally that if the 



scaling dimension of X2r+2 is positive then the asymp- 
totics in Eq. (Is]) gives a subleading correction to the reg- 
ular analytic part, Pex ~ const L'^ u^, coming from the 
high energy (ultra-violet) contribution to the susceptibil- 
ity. We will discuss its importance in the next section. 

Other possible generalizations of the scaling law Eq. (fTl) 
involve studies of defect production in systems where the 
dynamics describes the passage through quantum crit- 
ical lines. A concrete example of such a situation oc- 



curs in the transverse-field XY model (Divakaran et al. 



2008 Mukherjee et al. 2008 ) 



Here the quench takes one 
through a gapless line where the critical point occurs at 
the same momenta (fc = for the present case) at each 
point on the line. A detailed analysis shows that in such 
cases, for critical lines with z — u = I, the defect den- 
sity still obeys a univers al scaling law albeit with a dif 



ferent power: 



,1/3 



(Divakaran et al. 



2008 



Mondal 



et al.\ |2009[ [Mukherjee et aL[|2008[). The second si t uation 



involves the 2D Kitaev model (Sengupta et al. 2008), 



where a quench once again involves the passage through 
a gapless line with an energy gap vanishing for different 
momenta at different points on the line. It can be shown 
that in such a case, the defect density scales as ^/v for 
2D Kitaev model instead of the expected n ~ w behavior 
for 2D systems with z = v = 1 { Sengupta et al. . 2008 ). A 



generalization of these results for linear quenches through 
critical lines with arbitrary z and v has also been worked 



out (Mondal et al. 2008 2009 Mukherjee et al. 2008 



Sengupta et al. 2008). Many other situations involv- 



ing anisotropic phase transitions and quenching through 
multi-critical points were analyzed in literature leading to 
various deviations from the scaling ([4|) ( Bermudez et al. 



2009 Deng et al. 2008 Sen and Vishveshwara 2010 ) 



In order to detect experimentally the density of ex- 
citations generated by passing through a QCP, it is 
evident that one should distinguish between situations 



where such excitations are long-lived quasi-particles (as 
for nearly integrable systems) or decay after being cre- 
ated (as for non- integrable systems). In the first case, 
the presence of excitations above the ground state could 
for example be detected by measuring correlation func- 
tions long after the quench. This has been shown for a 
Quantum Ising chain linearly tuned through its quan- 
tum critical point (Cherng and Levitov 2006). The 



presence of defects with respect to the ferromagnetic 
ground state lead to exponentially decaying correlations 
of the order parameter superimposed, for slow enough 
quenches, to characteristic oscillations with period scal- 
ing with the quench velocity. This second feature is 



observed for abrupt quenches as well (Sengupta et al. 



2004) and is a consequence of the integrability of the 
model (Rossini et al. 20101. If in turn we consider a 



generic non-integrable system, it is necessary to express 
deviations from adiabaticity in terms of quantities, such 
as the excess energy or the entropy generated by passing 
through the QCP, which are not sensitive to the decay of 
quasi-particles, but can still be related to the density of 
excitations created close to the quantum critical point. 
Energy can be unambiguously determined both exper- 
imentally and numerically for both integrable and non- 
integrable systems and its scaling with the rate of quench 
can be used to differentiate between different nonadia- 



Kehrein 2010 



batic regimes (Eckstein and Kollar 2010 Moeckel and 



Polkovnikov and Gritsev 2008). The ex 



cess energy or equivalently heat (Q) (IPolkovnikov 2008b ) 



generated during the quench process is in general univer- 
sal if the process ends near the critical point. The scaling 
of Q is associated with the singularity of the susceptibility 



X2r+i at the critical point (De Grandi and Polkovnikov 



2010), which implies that for relevant perturbations in 



the thermodynamic limit 



|((i+2)i'/(zi/r+l) 



(10) 



Unless one considers cyclic processes, the drawback of 
using heat as a measure of non-adiabaticity is that it is 
hard to separate it from the adiabatic part of the energy 
change, corresponding to the limit u — >■ 0. Moreover, 
if the position of the QCP is not exactly known, the 
heat becomes sensitive to the nonuniversal details of the 
spectrum at the final point of the evolution. A way out 
could be to measure the higher moments of the energy or 
the whole distribution function of the energy, connected 



to the statistics of the work ( Silva 2008 Talkner et al. 



2007) in finite size systems. In particular, in the case of 



abrupt quenches close to critical points the statics of the 



work is characterized by sharp edge singularities ( Paraan 
and Silva 2009 Silva [2008 ). A related natural measure 
of non-adiabaticity is obtained by focusing on the entropy 
since entropy is conserved only for slow (adiabatic) pro- 
cesses, while is expected to increase as the system passes 
through the QCP. Moreover entropy production can be 
detected experimentally in certain systems, e.g. in cold 
atoms by driving the system to the weakly interacting 
regime, where the relation between entropy and energy 



is known (Luo et al 



2007). Theoretically, the quantifi- 



cation of entropy production in a closed quantum system 
is rather subtle. Indeed, the von Neumann's entropy of 
the entire system, being conserved throughout unitary 
evolution (Landau and Lifshitz^ ^1980j , cannot be a good 
characterization of deviations from adiabatic dynamics. 



However, the concept of diagonal-entropy (Polkovnikov 
[2008a ), defined as S'rf = - Y.n P^n ln(p„„), where /9„„ are 
the diagonal matrix elements of the density matrix in 
the instantaneous basis, avoids this difficulty. In station- 
ary systems, the diagonal entropy is nothing but the von 
Neumann's entropy of the time averaged density matrix, 
also called diagonal ensemble. It is clear that the diago- 
nal entropy is generated only due to nonadiabatic transi- 
tions and thus satisfies the key requirement of the ther- 
modynamic entropy: it is conserved for adiabatic pro- 
cesses, and can only increase or stay constant in closed 
systems if the initial state is stationary (iPolkovnikovl 
[2008a). For initial equilibrium states the diagonal en- 
tropy also satisfies fundamental thermodynamic relation: 
dE = TdS — FdX, where F = —{dxH) is the generalized 
force. For particular noninteracting models, the scaling 
of the diagonal entropy was found to be the same as 
that of the density of quasi-particles Q 



et al. 2010a|b Mukherjee et al. 



2008). I t is also possi- 
ble to analyze the entanglement entropy (Calabrese and 



(De Grandi 



Cardyl [20041 [2005] |Refael and Moore, .2004, .Vidal et al. 



If 



2003), i.e. the von Neumann's entropy of the reduced 



density matrix of a part of the system, and in particular 



at its time evolution following a quench (Cincio et al. 



2007 Pollmann et al. 2010 Sengupta and Sen 2009). 



For specific ID spin systems it was found that the entan- 
gleme nt entropy scales logarithmically with the quench 
time ( [Cincio et oH |2007[ [Pollmann e^"aH|2010[ ). Notice 
however that, at the moment, it is unclear how one can 
measure entanglement in many body systems and the en- 
tangle ment entropy in par ticular (se e some su ggestions in 
Refs. ( [Klich and Levitov[ [2009, Klich et al.\ |2006) ) and 
what its relation with the thermodynamic entropy is. 

Finally, another interesting question that has received 
attention is the connection between microscopic dynam- 
ics and thermodynamics in the semiclassical limit. In 
general, in classical systems there is no simple analogue to 
the instantaneous energy levels, the key concept in anal- 
ysis of quantum systems. Such analogue, however, does 
exist in the case of periodic motion. Then in the semiclas- 
sical limit the stationary levels are found from the Bohr 
quantization (or more accurately from the WKB approx- 
imation (Landau and Lifshitz 1981 1), which states that 



the reduced action in the stationary orbit should be quan- 
tized. In classical mechanics it is known that the reduced 
action is an adiabatic invariant, i.e. it is conserved for the 
adiabatic evolution (Landau and Lifshitz 1982) . From 
the previous discussion applied to quantum systems we 
can deduce that near singularities like second order phase 
transitions, conservation of adiabatic invariants should 



be violated and this is in deed the case (Landau and Lif- 
In Refs. ([Altland et al. |2009| [Itin and 



shitz 1982) 



Torma 2009a|b ) slow dynamics was analyzed for a par- 
ticular many-body generalization of the Landau-Zener 
model (closely related to the Dicke model) in the semi- 
classical limit. It was found that the nontrivial power 
law scaling of the number of excitations in this system 
(similar to Eq. ([I])) follows from the changes of adiabatic 
invariants near the singularity, which in turn corresponds 
to a quantum critical point in the thermodynamic limit. 
It is interesting that quantum fluctuations in this prob- 
lem entered only through the initial distribution of the 
adiabatic invariants but not through the equations of mo- 
tion. The corresponding truncated Wigner approxima- 
tion turned out to be very accurate in all regimes of the 



dynamics (Altland et al. 



2009 Kiegel 2009). It would 



be very important to understand precise connections be- 
tween transitions among microscopic energy levels in the 
quantum case and changes of suitable generalizations of 
adiabatic invariants in the classical limit. 



C. Slow dynamics in gapped and gapless systems. 

Up to now we have discussed the universal dynam- 
ics resulting from the variation of a control parameter 
A through a quantum critical point. However the dy- 
namics of interacting quantum systems has interesting 
regimes even when the system is fully gapped or gapless 
for the entire duration of the protocol. The classifica- 
tion of these regimes is important in order to understand 
dissipation and to develop optimal protocols minimizing 
non-adiabatic effects. Many of these questions are cur- 
rently a subject of intense theoretical research in different 
contexts, from quantum computation to transport. 

The general formulas Eq.([6|)-(l7|, which describe the 
density of excitations P^xiv) generated by a variation 
of the control parameter, tell us that if the system re- 
mains fully gapped throughout the evolution, then Pox 
and riox will have a quadratic scaling with v whenever 
the susceptibility X 2r+2 evaluated at the initial and final 
couplings is finite ( De Grandi and Polkovnikov[ 2010 



Rigolin et al. 2008 ) . A similar argument shows that the 
heat Q is also quadratic in v if X2r+i is finite. This 
quadratic scaling is characteristic of any quantum sys- 
tem. Let us point that in the standard Landau-Zener 
problem in the slow limit A^ ^ z; if we start in the 
ground state at t -> — cx) and let the system evolve up 
t — >■ -foo, the transitions to excited states are exponen- 
tially suppressed as a result of destructive interference 



between multiple transitions ( Vitanov 1999 Vitanov and 



Garraway 1996). At the same time in the intermediate 



stages of the evolution the transition probability reaches 
much higher values which scale only quadratically with 
the rate v. For example, if one considers a process which 
starts at to — >■ — oo the transition probability to the in- 
stantaneous excited state at the mome nt t in the slow 
limit u/A^ ^ 1 can be approximated by (Vitanov 1999 1: 



Pe 



v^A^ 



16(.g(i)2 + A2)3- 



(11) 



If i > there is an additional exponential term which 
leads to Eq. (Is]) at t — )■ +cx). If the process starts at 
to = exactly in the symmetric point, where g{0) = 
the quadratic asymptotics (111 is also recovered (Cuc- 



chietti et al. 20071. This scaling occurs as a result of 



the discontinuity of the first derivative of g{t) at the 
moment where the process starts or ends or following 
a discontinuity in any other point of the evolution (see 
e.g. Refs. ( [Damski and Zurek[ |2006[ [Divakaran et al.\ Falicov-Kimball model ^Eckstein and KollarJ 



and what sets the scaling of various quantities. 

A close inspection of the adiabatic susceptibility shows 
that in general the analytic (quadratic) part of the heat 
and energy of excitations on v comes from the high en- 
ergy (or ultra-violet) part of the spectrum, while the 
non-analytic part comes from low energies. This was in- 
deed shown to be the case in several situations, from 
the sine-Gordon model (De Grandi et al. |2010b), to the 



2010[ ) for particular cases). Likewise if there is a discon- the turning on interactions in a Fermi liquid ([Moeckel 



tinuity in the second order derivative of g{t) asymptoti- 
cally the transition probability in the LZ problem scales 
quadratically with acceleration. More generally for the 
protocol g{t) = go + v{t — t^Y /r\ 9{t — tp), where 9{t) 
is the step function, one can show that ( |De Grandi and 
Polkovnikovl 120101) 



Pcx(t ^ OO) 



2^2 



16(g2 + A2) 



2r+l ■ 



(12) 



As we discussed in the previous section this formula ap- 
plies even to sudden transitions (r = 0) where it reduces 
to the result of the ordinary perturbation theory. The 
same expression applies to the reverse process. If both 
the initial and final couplings are finite then the result- 
ing transition probability is asymptotically determined 
by the sum of probabilities associated with discontinu- 
ities of derivatives of g{t) at the initial and final times 
of the evolution plus additional interference terms which 
highly oscillate in the slow limit. Let us point out that 
in the LZ problem (and in general in gapped systems) 
one can suppress power law asymptotics of the transition 
probability by starting and ending the protocol infinitely 
smoothly, e.g. g{t) ~ go + 5i exp[-T/(i - to)]9{t - to). 
In this case only the non-analytic term in the transition 
probability survives and we are back to Eq. ([3]) where v 
is the time derivative of g{t) near the symmetric point 
where g{t) = 0. 

In gapless systems the situation becomes qualitatively 
different. In this case the adiabatic susceptibilities can 
diverge leading to non-analytic dependence of the corre- 
sponding quantities on w, as in the case of the crossing of 
a quantum critical point. For example, it is straightfor- 
ward to see that for marginal perturbations in a generic 
gapless phase the scaling dimension of the adiabatic sus- 
ceptibility dim[x2r+2] = d — 2zr. It becomes negative in 
low dimensions d < 2zr leading to a non-analytic scal- 
ing of the density of excitations with v. Thus depending 
on dimensionality in gapless systems one expects at least 
two different regimes of the response of the system to 
a slow external perturbation: analytic and non-analytic. 



These regimes were first suggested in Ref. (Polkovnikov 
and Gritsev 2008[ ) together with a third regime where 
adiabaticity is violated in the thermodynamic limit and 
Q or TT-cx become proportional to a power of the system 
size or some other large length scale associated with some 
irrelevant operator. In this regime, which can be realized 
in low-dimensional bosonic systems, the scaling Eq. (fl]) 
is violated. At the moment it is unclear how generic it is 



20101 and 



and Kehrein 2010 1 . As we pointed above the ultra-violet 



transitions can be suppressed by avoiding discontinuities 
in A(i) and its derivatives. However, this is not neces- 
sarily the case for the low energy non-analytic contribu- 
tion. To see this we need to reexpress the excess energy 
(or density of excitations) in terms of the total time of 
the process r. Doing this it was found that in an in- 
sulating, gapped phase, the details of the protocol are 
important and smoother protocols lead to a suppression 
of non-adiabatic effects, while in a gapless phase making 
A(i) smoother does not affect the heating (Eckstein and 
Kollar 2010 1 . This result can be again understood by an- 
alyzing the scaling dimension of the susceptibility X2r+i- 
According to our discussion for generic gapless systems 
its scaling dimension is negative when d + z < 2zr. 
Then Q ^ yi'i+^)/^^ (this result immediately follows from 
Eq. ( 10 1 by taking the limit v — )■ oo). On the other hand, 
V is related to the total quench time as u ^ 1/t^. Thus 
we see that in this case Q ~ l/r'^'i+^)/^^ i.e. indeed in- 
dependent on r. On the other hand for positive scaling 
dimension of X2r+i which is the case for d + z > 2zr 
and which is always true in gapped systems we have 
Q r^ v^ ~ l/(r)^''. Since in the adiabatic limit t is large 
we see that indeed the heat can be suppressed by increas- 
ing r and making A(t) smoother. An interesting open 
question is finding an optimal protocol for minimizing the 
non-adiabatic effects within given time r. It is plausible 
that the optimal power is determined by a vanishing scal- 
ing dimension of the corresponding adiabatic susceptibil- 
ity X2r+i- The questions of finding protocols minimizing 
non-adiabatic effects for gapped systems (with possibil- 
ity of crossing isolated quantum critical points) were also 
addressed by approximately minimizing the transition 
probability and identifying the Riemannian metric ten- 



sor underlying the adiabatic evolution ( Rezakhani et al. 



2010 2009) 



Studying the optimization of the protocol 
taking a system through a QCP it was found the optimal 
exponent of X{t) ^ \t\^sign{t) near the QCP scales log- 
arithmically with the quench time, r oc ln(r) (Barankov 
and Polkovnikov] 20081. This result was al so extended 
to systems with external confining potential ( [Collura and 
Karevskil 120101) . 



D. Effects of finite temperature. 

In the discussion above we always implicitly assumed 
that the system is initially prepared in the ground state. 



9 



An interesting and genuine question is how finite tem- 
perature effects modify the picture. In isolated systems 
temperature enters through initial conditions: the sys- 
tem is prepared in the initial finite temperature equilib- 
rium state and is then dynamically driven out of equi- 
librium. How is the response of the system affected 
by the initial thermal fluctuations ? One naturally ex- 
pects that while the transitions to high energy states 
(quadratic in v) will not be affected by small temper- 
atures in the system, the transitions to the low energy 
states, which determine the non-analytic contribution to 
heat and density of excitations, will be very sensitive to 

|2010a|b[ ) (see 



temperature. In Refs. ( De Grandi et al. 



also Ref. (Gritsev and Polkovnikov 2010)) studying a 



particular Sine-Gordon model in the two limits where it 
could be mapped to free bosons and free fermions, it was 
shown that the statistics of quasi-particles enters the scal- 
ing of both Q and riex making dynamics more adiabatic 
(compared to the zero temperature case) for fermions due 
to Pauli blocking and less adiabatic for bosons due to 
Bose enhancement. These results were not yet extended 
to generic interacting systems. 

Another aspect of thermalization, the influence of the 
coupling to an environment setting the temperature on 
the slow dynamics near quantum critical points, has been 



studied in Refs. (Fubini et al 



2007 



Patane et al.\ |2009a Patane et al. |2008[ [Patane et al. 
2009b). This setup allows one to analyze the effects of 



Mostame et al. 



2007 



thermal smearing and of dephasing/dissipation on the 
dynamics of a quantum critical system. Using a combi- 
nation of kinetic equations and scaling arguments it was 
found that in this situation the excess energy has two 
universal contributions, one still given by Eq. (10 1, while 
the second involving a universal power of temperature 
replacing the universal power of v ( Patane et al. 2008 ) . 



E. Open problems 

While the physics described above is definitely an im- 
portant example of the emergence of universality in the 
dynamics of interacting quantum systems, it is evidently 
a piece, albeit important, of the puzzle that has to be 
composed in order to understand to which extent the 
standard concepts of statistical physics can be applied to 
nonequilibrium problems. Understanding the meaning 
of relevance or irrelevance of a perturbation in generic 
nonequilibrium processes, extending the notion of uni- 
versality to nonequilibrium systems, as well as the con- 
cept of renormalization group, is a task (or dream) that 
certainly requires the solution of many specific problems, 
and a close comparison between experiments and theory. 

So far most of the theoretical research focused on ana- 
lyzing slow dynamics for global quenches, where the ex- 
ternal perturbation couples to the whole system. How 
these results can be extended to local or spatially nonuni- 
form perturbations is an open question. At one extreme 
limit, one can imagine performing a quench only locally. 



Then the rest of the system could be seen as a thermal 
bath. The analysis of a special case of dynamics of a 
transverse field Ising model where the tuning parame- 
ter linearly depends both on time and space has shown 
that excitations are generally suppressed by nonunifor- 



mity of the tuning parameter (Dziarmaga and Rams 
2010). This suggests that quantitative and qualitative 



differences may emerge when some of the symmetries of 
the system, e.g. translational, are broken in the quench 
process. 

Another important issue concerns the connections be- 
tween adiabaticity in thermodynamics and microscopic 
dynamics. One of the consequences of the thermody- 
namic adiabatic theorem is that no heat can be gen- 
erated in an isolated system during an infinitcsimally 
slow process. More generally according to the second 
law of thermodynamics in the Thompson's (Kelvin's) 
form for any cyclic process the system can only in- 
crease its energy, i.e. the heat should be always non- 
negative as long as one starts in equilibrium. This 
statement, which is obvious if the system is initially in 
the ground state, has been proven microscopically for a 
class of passive initial states ( Allahverdyan and Nieuwen- 
huizen[ [20021 [2005) [Boksenbojm et al.\ |2010| [Thirring 



2002), whose initial density matrix is stationary (diag- 
onal) and monotonically decreasing function of energy: 
(pn — Pm){^m — Cn) > 0. This Statement also directly 
follows from analyzing transitions between microscopic 
energy levels ( [Polkovnikov 2008b). Likewise many state- 
ments of thermodynamics related to behavior of entropy 
including the second law and fundamental thermody- 
namic relations are recovered using the concept of diago- 
nal entropy ( Polkovnikov , 2008a| ) . At the same time there 
are many open questions remaining: what are the time 
scales involved in the definition of adiabaticity ? How 
one can microscopically define adiabatic time scales in 
interacting systems and why these time scales are much 
shorter than inverse distance between many-body levels 
(see e.g. discussion in Ref. (Bahan 1991)) ? And finally, 
what is the role of integrability in nonequilibrium ther- 
modynamics? These questions are closely connected to 
the microscopic origin of conventional dissipation, which 
in turn is also very likely related to the combination of 
nonadiabatic creation of the elementary excitations and 
their following relaxation or dephasing. From the discus- 
sion above, we can anticipate anomalous dissipation near 
critical points and in gapless low-dimensional systems. 



III. EFFECTS OF INTEGRABILITY AND ITS 
BREAKING: ERGODICITY AND THERMALIZATION. 

Let us now turn to one of the most natural questions 
to be addressed when studying the dynamics of a closed 
many-body quantum system: are interactions within the 
system sufficient to make the system behave ergodically ? 
If we focus on local degrees of freedom, e.g. a few spins in 
a spin chain, can the rest of the system be always thought 
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as an effective thermal bath? And if this is not possible, 
are there some observable effects on the system dynam- 
ics? While these q uestions are definitely connected to 
quantum ergodicity (Goldstein et al. 2010), a topic with 



a long history dating back to the early days of quantum 



mechanics (Deutsch 1991 Mazur 1968 



von Neumann 



l929l|Paun and Fierz||1937[|Peres[[T984[ |Srednicki[ |1994 
Suzuki 1971), the past few years have brought a great 



deal of progress in the context of closed many-body sys- 
tems. The main motivation came from recent experi- 
ments on low dimensional cold atomic gases described 
in some detail in Sec. IV in this Colloquium (Greiner 



et al. 2002b Kinoshita et al. 2006). The experimental 



availability of essentially closed (on the time scales of ex- 
periments) strongly correlated systems together with the 
awareness of the conceptual importance of these issues 
in a number of areas (e.g. transport problems, many- 
body localization, integrable and non-integrable dynam- 
ics) have stimulated a lot of interest on quantum ther- 
malization. Below we will give a synthetic view on a 
number of recent important developments on this sub- 
ject, starting with the discussions of the general concepts 
of ergodicity and thermalization, and then moving to the 
discussion of many-body systems and integrability. 



A. Quantum and classical ergodicity. 

While the idea of ergodicity is well defined in classical 
mechanics, the concept of quantum ergodicity is some- 
what less precise and intuitive. Classically, an interacting 
system of N particles in d dimensions is described by a 
point X in a (2 d A'')-dimensional phase space. The intu- 
itive content of the word " ergodic" , i.e. the equivalence of 
phase space and time averages, can be then formalized by 
requiring that if we select an initial condition Xq having 
initial energy H{Xq) — E, where H is the Hamiltonian 
of the system, then 

6{X~X{t))= hm l^ f dtS{X-Xit)) = p,,,,{E), 

(13) 
where Pmc{E) is the microcanonical density of the system 
on the hyper-surface of the phase space of constant en- 
ergy E^ and X{t) is the phase space trajectory with initial 
condition X^. Of course if this condition is satisfied by all 
trajectories, then it is also true for every observable. We 
immediately see that in order to have ergodicity, the dy- 
namics cannot be arbitrary: the trajectories X{t) have 
to cover uniformly the energy hyper-surface for almost 
every initial condition Xq. 

The most obvious quantum generalization of this no- 
tion of ergodicity is arduous (von Neumann 1929). Let 



us first of all define a quantum microcanonical density 
matrix: given a Hamiltonian with eigenstates | ^Pq) of 
energies Ea, a viable definition of the microcanonical 
ensemble is obtained by coarse graining the spectrum 
on energy shells of width 5E, sufficiently big to con- 
tain many states but small on macroscopic scales. De- 



noting by T-L{E) the set of eigenstates of H having en- 
ergies between E and E + SE, we define p„ic{E) = 
J2a<£H(E) 1/-^ I *Q)(^a I; whcre Af IS the total num- 
ber of states in the micro-canonical shell. Let us now 
ask the most obvious question: given a generic initial 
condition made out of states in a microcanonical shell, 

I *o) = E 



ae«(B) 



I *c<), is the long time average 
of the density matrix of the system given by the micro- 
canonical density matrix? The answer to this question 
for a quantum system is, unlike in the classical case, al- 
most always no, as J. von Neumann realized already in 
1929 (von Neumann 1929). More precisely, if we assume 



the eigenstates of the system not to be degenerate, the 
time average is 



*(i))(*(i)|-5]|c„n*a)(*«h 



Pdia 



(14) 



where | \l/(t)) is the time evolved of | V^o). This object is 
known in the modern literature as the diagonal ensem- 
ble ( [Rijoll [2009| |Rigol et al.\ [20081 [2007| ). Notice now 
that the most obvious definition of ergodicity, i.e. the 
requirement pmc = Pdiag, implies that | Ca P= 1/Af for 
every a, a condition that can be satisfied only for a very 
special class of states. Quantum ergodicity in the strict 
sense above is therefore almost never realizable (Gold- 



stein ei aL 2010 von Neumann 1929). 



Our common sense and expectations, which very fre- 
quently fail miserably in the quantum realm, make us 
nevertheless believe that, in contrast with the arguments 
above, macroscopic many-body systems should behave 
ergodically almost always, unless some very special con- 
ditions are met (e.g. integrability). The key to under- 
stand ergodicity in therefore to look at quantum systems 
in a different way, shifting the focus on observables rather 



than on the states themselves (Mazur 1968 von Neu 



mann 



1929 Peres 



1984 ). Given a set of m,acroscopic 
observables {M^Ya. natural expectation from an ergodic 
system would be for every | ^o) on a microcanonical shell 

nE) 

(*(i) I Afyit) I *(t)) -^t^+oo Tr[M^/5^e] = (M^),„e(15) 

i.e. that looking at macroscopic observables long after the 
time evolution started makes the system appear ergodic 
for every initial condition we may choose in 'H(E). One 
needs a certain care in defining the infinite time limit 
here, since literally speaking it does not exist in finite 
systems because of quantum revivals. A proper way to 
understand this limit is to require that Eq. ( 15 ) holds in 



the long time limit at almost all times. Mathematically 
this means that the mean square difference between the 



LHS and RHS of Eq. (15) average d over long tim es is 
vanishingly small for large systems dReimann 2008 ) . To 



avoid dealing with these issues ergodicity can 
using the time average, i.e. requiring that 



3e defined 



(*(t) I M^{t) I *(i)) = Tr[Af^/5diag] = (Mp)^,. (16) 

Notice that if the expectation value of Mp relaxes to a 
well defined state in the sense described above, this state 
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will coincide with the time averaged state and the two 
definitions of ergodicity Eq.(15)-(16) will be equivalent. 



If the conditions above are satisfied then in loose terms 
Prac can be considered as equivalent to Pdiag • J- von 
Neumann proved that if the system satisfies some very 
natural requirements (e.g. absence of resonances), and 
the set {Mp} is constructed in such a way as to define 
macrostates of the system, which obviously requires the 
observables to be coarse grained on the various micro- 
canonical shells 'H(-E) and mutually commuting, then a 
form of ergodicity is observed (sometimes referred to a 
normal typicality). In particular, for every | ^o) and 
almost every set {M^} the diagonal and microcanoni- 
cal ensembles ar e equivalent ( Goldstein et al. . 2010[ |von 



Neumann 1929). More recently it was proven that the 



whole density matrix of a small subsystem of a bigger 
system which is placed initially in a typical eigenstate 



is described by the canonical ensemble (Popescu et al. 
[2006) . In Ref . (Gogolin et al., 2011J these results were 
further extended to the problem of measurement and de- 
coherence. Particular care is nevertheless needed in relat- 
ing these statements to the dynamics and thermalization 
of actual many body systems, since physical initial con- 
ditions in quenched system almost never correspond to 
eigenstates of a new Hamiltonian. 



B. Nonergodic behavior in integrable systems: the 
generalized Gibbs ensemble 

While the statements above are very general, their ap- 
plication to specific systems is not at all straightforward. 
Looking at a concrete many-body system, it is of primary 
interest not just to find out whether in principle a set of 
macroscopic observables that behave ergodically exists, 
but whether specific and natural observables, such as the 
magnetization for spin chains, density for cold atomic 
gases, or various correlation functions behave ergodically 
or not. In this respect, exp eriments tell us that er godic- 
ity is not at all guaranteed ( Kinoshita et al. 2006 1 if the 
closed system is integrable or nearly integrable. While 



this fact was expected pBarouch etal. 1970 Girardeau 



1969 1970 Mazur 19681 Suzuki 1971 1 recent research on 



the dynamics of integrable systems has focused on find- 
ing ways to predict the asymptotical states taking into 
account integrability, i.e. the presence of many constants 
of motion. 

Let us discuss how this can be done qualitatively us- 
ing the simplest example of integrable system, a periodic 
harmonic chain of finite length described by the Hamil- 
tonian 



H 



A/-1 

E 



2m 



-{Xj -Xj + iY' 



(17) 



where Xj are deviations of particles from the equilibrium 
positions and pj are their momenta; we use the identi- 
fication xm = xq. Let us imagine that initially we de- 
form the system in a particular way and ask how this 



deformation evolves with time. We note that since this 
is a harmonic system described by linear equations of 
motion the following analysis also applies to quantum 
systems. From elementary physics we know that the ini- 
tial deformation splits into normal modes characterized 
by the quasi-momenta g„ — 2T:n/M, where n is integer 
n e [0, M — 1], and the dispersion ujq = 2i/smq/2. This 
system obviously does not thermalize even at long times 
because there is no energy exchange between modes. 
This does not imply though that it can not reach a well 
defined asymptotic state in the long time limit (Barthel 
and Schollwock[ |2008| [Cramer and Eisert[ |2010[ ). To il- 
lustrate this point consider for example the displacement 
of the j-th atom at time t after some initial displacement: 



Xj{t) 



1 

/M 



M-l 



E^«"We^'" 



(18) 



Xq{t) = Aq cos[a;^t], where Aq is a complex amplitude de- 
termined by initial conditions (for simplicity we assumed 
an initial stationary state). Let us now analyze quali- 
tatively the dynamics of this system. At short times, 
provided that the initial modulation is smooth and only 
modes with small momenta [q <IC tt) are excited we can 
linearize the spectrum ujq « uq. Clearly in this case we 
recover periodic motion of the wave-packet with a period 
equal to the ratio of the system size and sound velocity: 
T = M/v. This persistent motion is characteristic of the 
absence of any relaxation. However, as time gets longer 
deviations of the dispersion from linear become more im- 
portant. In particular, when t*[(jjn^i -\-ijJn-i — 2ujn) ~ 1, 
where ujn is the central frequency of the wavepacket, cor- 
relations between phases among the different modes are 
lost and they can be treated as essentially random num- 
bers. For our model t* ~ M"^ /uJn. At long times t ^ t* 
the different momentum modes become uncorrelated and 
the system reaches the asymptotic stationary state in a 
sense we defined earlier (it can be found close to that 
state at almost all times). For this asymptotic state the 
only relevant information about the initial conditions is 
encoded in the M mode amplitudes \Aq\ or equivalently 
in their squares \Aq\'^ proportional to the occupancies of 
the modes of energy Eg, which are the integrals of mo- 
tion. Note though that there are special modes corre- 
sponding to momenta q and —q which are exactly degen- 
erate. The correlations between Aq and A-q = A* thus 
never disappear and in general one needs to fix M addi- 
tional constraints representing the relative phases of the 
complex amplitudes. For example, if the initial configu- 



ration IS symmetric Xj 



then Aq is real, meaning 



that the phases of all modes are identical. Then it is 
easy to check that with this constraint (a;?(t)) acquires 
spatial dependence on j even in the long time limit. This 
dependence can not be recovered by fixing only mode 
occupancies. Only when these phases are unimportant 
e.g. they average to zero or if there are no degenera- 
cies between the normal modes the asymptotic state is 
fully fixed by the integrals of motion. Thus in contrast 
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with ergodic systems where only the energy needs to be 
fixed, the long time behavior of our integrable model can 
be reproduced by fixing the M integrals of motion and 
possibly ^ M other constraints if there are degenera- 
cies. While the number of commuting (local) integrals of 
motion is large, equal to the system size M, it is vastly 
smaller than the total number of states which scales with 
M exponentially. 

Let us now see how these considerations are transposed 
in many-body systems, focusing on another simple inte- 



grable model, the Quantum Ising chain (Sachdev 1999) 
described by the Hamiltonian: Hq = — Yli'^i'^i+i+9'^i- 
Here af '^ are the spin operators at site i and g is the 
strength of the transverse field. This model gives one of 
the simplest examples of quantum phase transition, with 
a quantum critical point at ^c = 1 separating two mutu- 
ally dual gapped phases, a quantum paramagnet {g > gc) 
and a ferromagnet (g < gc). 

In the Quantum Ising chain the local transverse magne- 
tization, M'^ = ^,- gf is a non-ergodic operat or ( Barouch 



et a. 1970 



Girardeau 1969 1970 Mazur 



1968). To 



see this, it is useful to employ a Jordan- Wigner trans- 
formation that reduces the problem to a free fermion 
model ( [Sachdev 19991. In terms of the fermionic op- 
erators Cfc relative to modes of momentum k = 2TTn/L 
the Hamiltonian takes the form 

H = 2^(g-cos(fc))(4cfe-c_fcC^_,) 



k>0 



isin(fc)(c[clj, -c_fcCfc), 



(19) 



Under this mapping the transverse magnetization be- 
comes M^ — ~'2J2k>o(^k^k — c-fcclj,). The eigenmodes 
7fc of energy Ek — '2\/{g ~ cos{k))^ + sm{k)^ diagonal- 
izing the Hamiltonian are related to the fermionic op- 
erators Ck by a Bogoliubov rotation, Ck — cos{Ok)jk — 
isin{0k)T_i., with tan(26'fe) = sin(fc)/((7 — cos(fc)). In 
Heisenberg representation the operators jk acquire sim- 
ple time dependence: jk{t) — 7fe(0) exp[— iiJfct]. As in 
the previous problem of a harmonic chain if the energies 
Ek are incommensurate, at sufficiently long times differ- 
ent momentum modes become statistically independent 
from each other. This statement does not apply to modes 
with opposite momenta k and —k which have identical 
energies. However, if these correlations are not important 
then in the long time limit (see below) each mode can be 
characterized by the conserved quantity ri/c = {'Tf.lk)- Let 
us now consider unitary dynamics of the transverse mag- 
netization starting with a generic initial condition | ^o)- 
The time evolution of the operator Af^(i) expressed in 
terms of the eigenmodes of the Hamiltonian is 



M-(t) = -2 Vcos(20,)(7^7, -7_,.7l,) 



fc>0 



+ i sin(20fc)(7_fc7fce-2^^'=* - 7T7 



t^,t ^2iEkt\ 



(20) 



In the long time limit only the diagonal terms in the sum 
survive, while the off diagonal ones, describing creation or 



destruction of two fermions average to zero. Therefore for 
any initial condition | ^g) the asymptotic value attained 
by the transverse magnetization is 



(M-(t)) 



-2^cos(20,)((7i.7fc) - {i-kllk))- (21) 

fe>0 



This asymptotic value is therefore perfectly described by 
the set of the occupation numbers n^. 

The result above leads one to conjecture that the 
asymptotic state is described by a Gibbs-like statistical 



ensemble of the type ( Rigol et al. 2007 1 



e^T^k^kikik 



PG 



(22) 



where the Lagrange multipliers Afc are fixed by requiring 

that Uk = (*o I 7I7/C I *o) = Tr[pG7l7/c] = (7l7fe)G- 
The ensemble defined above in Eq. ([22]) can be seen as a 
particular case of the ensemble 



E„A„7„ 



Pg 



(23) 



known as the generalized Gibbs ensemble (GGE) or 
the maximum entropy ensemble and introduced by 



Jaynes ( Jaynes 1957 ) to describe the equilibrium state of 
a system possessing N constant of motions Ik- A recent 
conjecture ( ,Rigol et al.^ ^2007) proposed to use the GGE 
to describe the asymptotic state of a generic quantum 
integrable model. This proposal had however to face two 
obvious subtleties. It is first of all to be specified how 
to choose the Ik in Eq.(23). Indeed, if all constants of 



motion would be admissible, including non-local ones, 
then one would obviously and tautologically describe the 
asymptotic state (both for integrable and non-integrable 
systems), as one can easily see by choosing as Ik the pro- 
jectors onto the eigenstates of the Hamiltonian. The way 
out comes however by observing that in standard ther- 
modynamics the Gibbs ensemble emerges for small sub- 
systems from the assumption of statistical independence 
between sufficiently big subsystems. In this derivation 
the additivity of a conserved quantity - energy - plays a 
crucial role. This is the reason why the probability of 
a given configuration is exponential in energy and not 
e.g. in energy squared (Kardar 2007). Similar argu- 
ments apply to any additive integrals of motion so that 
statistical independence and invariance of the ensemble 
to the choice of a subsystem of an integrable system puts 
strong constraints on the choice of the integrals of mo- 
tion in GGE when the latter is applied to subsystems 
of an integrable system. In this respect the average oc- 
cupation numbers of different momentum modes used in 
Eq. (22 1 become approximately additive for small subsys- 



tems. This approximate additivity of integrals of motion 
was recently discovered in Ref . ( Cassidy et al. , |2011 ) for 
another integrable system of one-dimensional hard-core 
bosons. In particular, it was noticed that the integrals of 
motion /„ and the lagrange multipliers A^ in that case 
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can be written as a smooth functions of a/N implying 
that in large systems the argument of the exponent en- 
tering Eq. ([23]) can be written in the extensive (additive) 

form: X^a ■^a-^a ~ i/o d^Ji.OKO^ where ^ = a/L. 
This suggests that GGE can be defined through a smooth 
function A(^), which replaces the temperature in the er- 
godic systems. 

There is a second subtlety in applying the GGE to 
quantum systems. Here the most natural definition of 
integrability is based on requirement that the system 
has well defined quasi-particles that maintain their iden- 



tity upon scattering (see Ref. (Gaux and Mossel 2011) 



for a more detailed discussion), i.e. scattering is purely 
elastic and there is no production of particles or dis- 



sipat ion associated to it (Mussardo 2009 Sutherland 
2004). This notion can be made precise in continuum 



integrable models, such as the Luttinger liquid or the 
Sinh-Gordon model, which can emerge as low energy de- 
scriptions of other integrable models, such as the critical 
XXZ chain and the Lieb-Liniger gas. In these systems 
it is natural to associate the I^ to the occupation num- 
bers of the quasi-particle states. More specifically, con- 
sidering a generic one dimensional relativistically invari- 
ant integrable system with say a single species of quasi- 
particles of mass TO, energy E = m cosh(0) and momen- 
tum p = TOsinh(0) {9 is the rapidity), the quasi-particles 
can be described by annihilation operators A{9) satis- 
fying the algebra A{e^)A{ej) = SiOi - 0j)A{ej)A{ei), 
where S is the S-matrix of the two particle scattering. 
Similar relations are valid for the products of creation, 
and creation- annihilation operators (Mussardo 2009). 



Since the Hamiltonian is by definition diagonal in A{0), 
H = J d0E{9)A\9)A{9), and every eigenstate can be 

written as | 6'i, . . . ,6i„) = At (6*1) At(6i„) | 0), with 

61 > ■ ■ ■ > 6n, in this case it is rather natural to postulate 
the form 



PG 



' f dex{e)A'*(e)A(9) 
Z ' 



(24) 



for the generalized Gibbs ensemble (Fioretto and Mus- 



sardo 2010). This ensemble is a direct generalization of 
the GGE for the Quantum Ising Model, where 5' = — 1. 
For this general class of integrable systems and a spe- 
cific class of translationally invariant initial states it was 
indeed shown the long-time limit of the average of local 
operators is well described by this ensemble ( Fioretto and 



Mussardo 2010). Such initial states can be written as 



I ^0) ^Ue-i'^^'^^''^'^^^^^^^'^~'^\ 



(25) 



which in turn are similar to the so-called integrable 



boundary states in statistical field theory (Ghoshal and 
[Zamolodchikovj 1994). Such states naturally emerge 
in experimentally relevant systems, for example when 
studying dephasing in split quasi- Id condensates (Grit 



sev et al.. 2007) or in the Quantum Ising model, when 



studying a quantum quench from a transverse field 7,^ to 
a transverse field 7/ ( Silva 2008 ) . 



A very interesti ng idea related to the GGE was sug- 
gested by Gurarie (Gurarie 19951 to explain the steady 



state of a driven nearly integrable system. It was shown 
that the steady state distribution of the wave amplitudes 
corresponding to different momenta (see Ref. (Zakharov 
et al. . 1992 ) for details) can be obtained by taking the 
probability density p ex exp[— _F], where i^ is a (complex) 
combination of the approximate integrals of motion found 
perturbatively. In terms of this ensemble one recovers the 
correct power law distribution of the amplitudes of waves 
with the momentum and other observables. 

Another view towards elucidating the validity of the 
generalized Gibbs ensemble has been pursued for special 



quenches in a ID Bose-Hubbard model (Cramer et al. 
20081) and in integrable systems with free quasiparti- 



cleslBarthel and SchoUwock 2008). It was shown that. 



upon tracing all degrees of freedom of the system outside 
a small region of space and under specific conditions, the 
local density matrix tends asymptotically to po- More 
recently a series of recent theoretical ( Flesch et al. 2008 ) 
and experimental (Trotzky et al. 20111 works on the 



dynamics of Bose-Hubbard models has proven the re- 
laxation of local observables in these system to a max- 
imum entropy ensemble consistent with the constraints 
of the dynamics. A hint towards the generalization of 
Pg for Bethe Ansatz integrable systems was proposed 
in Ref. (IBarthel and Schollwockl [2008|. The GGE was 



also tested in a number of mode 



s, from Luttinger liq- 



uids (Cazalilla 2006 lucci and Cazalilla' '2009) and free 



bosonic theories ( Calabrese and Gardy^ .2007), to inte- 



grable hard-core boson models ( Rigol et al. 20071 and 



Hubbard-like models (Eckstein and KoUar 



lar and Eckstein 



2008 Kol- 



2008 1. In all cases, it was shown to 



correctly predict the asymptotic momentum distribution 
functions for a variety of systems and quantum quenches. 
At this point is should be stressed that as we discussed 
before the GGE does not always give complete descrip- 
tion of the asymptotic state of the system. In the simple 
example of the harmonic chain we saw that for generic 
initial conditions it is necessary to specify 2A'^ real con- 
stants or N complex amplitudes in order to correctly de- 
scribe the asymptotic state even if we focus exclusively 
on local observables. For a quantum Ising chain, more- 
over, pq can be interpreted as a grand-canonical dis- 
tribution with an energy dependent chemical potential 
^Ji■k — Ek — ^k- It is evident now that if we consider 
the correlations of Suk — jljk — {iklk), the occupation 
numbers of different eigenmodes, the GGE necessarily 
predicts {Sn^Srik') = 0. Likewise the GGE predicts the 
correlators of the type {'^\{t)^-kit)) are always equal to 
zero. For a generic initial state | ^0) both statements are 
not necessarily true: by e.g. breaking translational in- 
variance in the initial state one could have (^0 I 5ni.5nk> \ 
*o) 7^ and (*o | lll-k \ *o) ^ 0. Notice that the mere 
survival of off-diagonal correlations of this type when the 
evolution starts with a non-translationally invariant state 
signals in a sense the integrability of the model, i.e. the 
existence of well defined quasiparticles 7^. Indeed, follow- 
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ing the argument of Gangardt and Pustilnik ( Gangardt 



and Pustilnik 2008 1 , if the Hamiltonian of the system is 



translationahy invariant but integrabiUty is broken, the 
off-diagonal correlators are expected to decay to zero for 
any initial condition, thereby restoring the translational 
invariance in the asymptotic state. Finally, notice that 
off-diagonal correlations might influence the asymptotics 
of physically relevant observables: a simple example is 
the asymptotic value of ((M^(i))^), which for a generic 
non-translationally invariant condition | ^I'o) cannot be 
predicted using the GGE. 

A very important open question is to understand un- 
der which general circumstances the GGE can be ap- 
plied. For free fermionic and bosonic syste ms the GGE 
was argued to hold for local observables (Barthel and 



Schollwock 2008 Cramer ei a/. 2008). For more general 



integrable systems this is not evident at all. For example 
in the case of the Quantum Ising chain the two signif- 
icant observables, the transverse magnetization af and 
the order parameter ct|, are local in the spin representa- 
tion. However, this locality does not translate directly to 
their representation in terms of the quasiparticles of the 
model: while af retains a local character in terms of 7i, 
erf does not. Will the asymptotic dynamics of any local 
operator be represented by the GGE, or just that of local 
operators in the quasi-particle fields ? Do the symmetries 
of the initial state play any role in this ? Answering these 
questions appears to be crucial to understand the role of 
integrability in the dynamics of many-body systems. 

Another very important question is whether all natu- 
ral observables of an integrable system behave necessar- 
ily nonthermally. The answer to it appears to be no, as 



pointed out recently ( Rossini et al. 2009). The key to un- 
derstand this issue seems to be again locality with respect 
to the quasi-particles diagonalizing the model. Thus in 
the Quantum Ising Model it was shown that while the 
transverse magentization is non-ergodic, the correlators 
of the order parameter az following a quench of the trans- 
verse field relax as in a thermal state with an effective 
temperature Tcff set by the initial energy of the system 
E = (*o(ft) I H{gf) I ^o(5*)>- At low Toff this relaxation 
appears to be universal, i.e. determined only by the low 
energy scattering properties of quasi-particles (Rossini 



et al. 2009). Analogous studies for an XXZ chain hint 



towards a different behavior of local and non-local opera- 
tors with respect to quasi-particles (Canovi et a?. |2011 ). 
The situation is much less clear for quenches with high 
effective temperatures, where the universal character of 
the low energy theory is lost (Barmettler et al. 2009). 



C. Breaking integrability: eigenstate thermalization. 

When integrability is explicitly broken with a strong 
enough perturbation one naturally expects ergo die be- 



2007 



havior to emerge for all observ ables ([Kollath et al. 
Maninana eTd] [2007{ |RigoH [2009| |Rigol et al 
Roux 2009 2010). The quest for the necessary con 



2008 



ditions for thermalization to occur (i.e. how strongly 
should integrability be broken, which spectral properties 
should the system display) is an important problem in 
many different fields, from mathematical and statistical 



physics to quantum chaos (Deutsch 1991 Peres 1984 



Rigoleia/. 2008 Srednicki 1994 1999). In classical sys- 



tems the intense research on this subject was stimulated 
by the study of dynamics of a nonlinear chain of cou- 
pled oscillators by Fermi Pasta and Ulam (FPU) ( [Fermi 
et al. May 1955 ) , where instead of thermalization regular 



quasi-periodic oscillations were observed. Later it was re- 
alized that the FPU problem is nearly integrable and that 
there is a finite threshold for the chaotic behavior ( Camp- 



bell et al. 2005 ) . In quantum systems the situation was 



far less clear: while different views on this issue emerged 
from time to time, the key towards a clear understanding 
of quantum thermalization appears to be linked to the 
emergence of quantum chaotic behavior ( [Peres 1984). 
In particular, it has been proposed that the emergence 
of thermal behavior is linked to the pseudo-random form 
of natural observables once represented in the eigenbasis 
of the Hamiltonian (Peres 1984). This observation has 



been made more precise by conjecturing that thermal- 
ization in quantum chaotic systems occurs eigenstate-by- 
eigenstate, i.e. the expectation value of a natural observ- 
able ($Q I A I ^c) on an eigenstate | ^q.) is a smooth 
function of its energy E^ being essentially constant on 
each microcanonical energy shell (Deutsch 1991 Sred- 



nicki 1994 1999). If this happens, then ergodicity and 



thermalization in the asymptotic state follow for every 
initial condition sufficiently narrow in energy (e.g. local- 
ized in a microcanonical shell), as one can easily under- 
stand using the diagonal ensemble. This hypothesis is 
known as eigenstate thermalization (ETH). 

In order to understand how eigenstate thermalization 
can emerge, let us consider a quantum gas of N par- 
ticles of mass m with hard-core interactions (Srednicki] 
1994). Srednicki pointed out that in the time evolution of 
this system starting with an initial condition | ^I'o) suf- 
ficiently narrow in energy, the momentum distribution 
will always relax to the Maxwell-Boltzmann distribution 
/mb(p) as long as the eigenstates of the system | '^a) can 
be considered as pseudo-random superpositions of plane 
waves, i.e. have a diffusive nature in phase space. This re- 
quirement should be satisfied as a result of the chaoticity 
of the system, the so-called Berry's conjecture. Calling 
X = (xi,X2, , . . . ,XAr) the coordinates of the particles 
and P = (pi, P2, • ■ • Vn) their momenta. Berry's conjec- 
ture states that the eigenstates have the form 



*„(X) =N I rfPA„(P)(5(p2 - 2mEa)e' 



p X 



(26) 



with Aq,(P) being pseudo-random variables with gaus- 
sian statistics, (A„(P)A/3(P')) = 5i5^'^^\V+I")/5{V^- 
P' ). Notice that we are disregarding the symmetrization 



of the wave function, see (Srednicki, 1994) for a discus- 
sion of this aspect. If the properties above are assumed 
it is easy to prove that on average in the thermodynamic 
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limit the momentum distribution function is 

(/(p)> = Idp2dp3 ...(|^„(p,P2,...)l') 

fMBip), (27) 



e -ImkT 



(27rTOfcT)3/2 

where the temperature is set by the equipartition law as 
Ea = 3/2NkT. Notice that this is expected to happen 
for every eigenstate of energy close to E^ , as required by 
the ETH. Hence thermal behavior will follow for every 
initial condition sufficiently narrow in energy. 

For generic many-body systems, such as Hubbard-like 
models and spin chains, the close relation between break- 
ing of integrability and quantum chaotic behavior is a 
known fact (Poilblanc et al. 19931. In particular, finite 



size many-body integrable systems are characterized by 
the Poisson spectral statistics while the gradual breaking 
of integrability by a perturbation leads to a crossover to 
the Wigner-Dyson statistics. The latter is typically asso- 
ciated, in mesoscopic systems or billiards, with diffusive 
behavior and can be taken as a signature of quantum 
chaos (Imry 1997). In many-body disordered systems 



the emergence of the Wigner-Dyson statistics was ar- 
gued to be an indicator of the transition between metallic 



(ergodic) and insulating (non-ergod i c) ph ases (Mukerjee 
et al.\ 2006 Oganesyan and Huse 2007). Inspired by 



these close analogies, recent studies gave a boost to our 
understanding of the crossover from non-ergodic to ther- 
mal behavior as integrability is gradually broken and of 
the origin of ergodicity/thermalization in systems sufh- 

Kollath 



ciently far fr om integrability ([Biroli et al. 
eTaTl [20071 [Manmana et al.\ |2007[ |Rigo 



2010 
2009 



> t 



Rigol 



et al. 2008 ). In particular, a careful study of the asymp- 



totics of density-density correlators and momentum dis- 
tribution function for hard-core bosons in Id showed that 
the transition from non-thermal to thermal behavior in 
finite size systems takes the form of a crossover con- 
trolled by the strength of the integrability breaking per- 
turbation and the system size (Rigol 2009). Moreover 



there is a universality in state to state fluctuations of 



simple observables in this crossover regime (Neuenhahn 



and Marquardt 2010), which goes hand-by- hand with 



less attention. In a series of studies of relaxation in 
fermionic Hubbard models subject to quenches in the 
interactions strength it has been argued that for suffi- 
ciently rapid quenches relaxation towards thermal equi- 
librium occurs through a pre-thermalized phase ( Moeckel 



and Kehrein 2008 20101. Similar two step dynamics oc- 



curs in quenches of coupled superfluids where initial fast 
"light cone" dynamics leads to a pre-thermalized steady 
state, which then slowly decays to the thermal equilib- 
rium through the vortex-antivortex unbinding ( |Math ey 



and Polkovnikov 20101. In Ref. (Burkov et al. 12007) a 



very unusual sub-exponential in time decay of correlation 
functi ons was predicted and later observed experimen- 
tally ( Hofferberth et al. 2007 ) for relaxational dynamics 



of decoupled Id bosonic systems. 



D. Outlook and open problems: quantum KAM threshold 
as a many-body delocalization transition ? 

The arguments above clearly pointed to the connection 
between thermalization in strongly correlated systems 
and in chaotic billiards. This analogy however , rather 
than being the end of a quest, opens an entire new kind 
of questions, which are a current focus of both theoreti- 
cal and experimental research. In particular, we do know 
that in a number of models of strongly correlated parti- 
cles eigenstate thermalization is at the root of thermal 
behavior (Rigol et al.\ 2008). What is the cause of eigen- 



state thermalization in a generic many-body system, i.e. 
the analogue of the diffusive eigenstates in phase space 
of Berry's conjecture ? And most importantly, while in a 
finite size system the transition from non-ergodic to er- 
godic behavior takes the form of a crossover, what hap- 
pens in the thermodynamic limit ? Is the transition from 
ergodic to non-ergodic behavior still a crossover or it is 
sharp (a quantum, KAM threshold ) ? 

At present, research on these questions has just 
started. An interesting idea that has recently emerged is 
that the study of the transition from integrability to non- 
integrability in quantum many-body systems is deeply 
connected to another important problem at the frontier 
of condensed matter physics: the concept of many-body 



an analogous transition from Poisson to Wigner-Dyson locahzation (Altshuler et al. 1997 Basko et al. 20061, 



level st atistics (Rigol and Santos 2010 Santos and Rigol which extends the original work of Anderson on single- 
2010a). When integrability is broken by sufficiently particle localization (Anderson 1958). We note that re- 



strong perturbation ergodic behavior emerges (Neuen- 
hahn and Marquardt[ [JUTOl [RigoU [20091 |Rigol and San- 



tos 



lated ideas were put forward in studying energy transfer 
in interacting harmonic systems in the context of large or- 



2010), which in turn appears to be related to the 



validity of the ETH ( [Rigol et al\ [2008[ ). In this con 



text, the anomalous, non-ergodic behavior of integrable 
models has been reinterpreted as originating from wide 
fluctuations of the expectation value of natural obscrv- 



ganic mo l ecules ( [Leitner and Wolynes[ [1996| [Logan and 
Wolynes 19901. More specifically, it has been noticed 



that a transition from localized to delocalized states ei- 



ables around the microcanonical average (Biroli et al. 
20T0|. 



ther in real space ( Pal and Huse 2010 1 or mo re gener- 
ally in quasi-particle space (Canovi et al. 2011 ) is closely 



All these statements apply to the asymptotic (or time 
averaged) state. So far the relaxation in time, in par- 
ticular in the thermodynamic limit, has received much 



connected to a corresponding transition from thermal to 
non-thermal behavior in the asymptotics of significant 
observables. For weakly perturbed integrable models, 
the main characteristic of the observables to display such 
transition appears to be again their locality with respect 
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to the quasiparticles (Canovi et al. 2011). This con- 



nection with many-body locahzation becomes more clear 
on the basis of a recently proposed way to quantify the 
transition from non-ergodic to ergodic behavior in many 
body systems (Olshanii and Yurovsky^ ^2009) . The au- 
thors consider an integrable Hamiltonian Hq with a weak 
non-integrable perturbation XV . Formulating essentially 
a generalization of the Berry's conjecture and making 
some additional assumptions they showed that the devi- 
ations from thermal behavior in the expectation value of 
observables can be quantified according to the formula 



{^{t) I A I ^{t)) - (A),„, « r;((vI/o(t) I A I vI/o(t))o - (^)mc), 

where | ^o(i)) = exp[— iHo*] I *o) is the time evolved 
state with respect to the integrable Hamiltonian, while 
I *(i)) = exp[-i(-Ho + AV)i] | *o) evolves with the 
non-integrable one. The key ingredient in this formula 
is the parameter 77, defined as the average over the 
microcanonical shell of the inverse participation ratio 
"^OL = Yl,n I (V'n I *q) I^, where | (^„) are the eigen- 
states of the integrable Hamiltonian 'Ho and | ^I'q) are 
the eigenstates of 'Hq -I- AV. Notice that when the system 
is close to integrability 77 ~ 1 but as the strength of V 
increases, r\ is roughly proportional to the inverse of the 
number of states M hybridized by the perturbation. 

Using this formula it is now possible to understand 
how many-body localization enters the scenario (Canovi| 
et al. 2011 Pal and Huse[ 2010): an abrupt transition 
at a certain Ac from localized to delocalized states in 
quasi-particle space corresponds to a sharp decrease of 
the inverse participation ratio rj from a value 0(1) to a 
value negligibly small and tending to zero in the ther- 
modynamic limit, essentially 0(1/A/'(A)), where A/'(A) is 
the total number of states in an energy window of width 
of the order of the matrix elements of the perturbation. 
This would lead to an abrupt transition from non-thermal 
to thermal behavior at Ac, a transition qualitatively cor- 
responding to the physics of a quantum KAM thresh- 
old. Notice that on the delocalized side of this transi- 
tion the eigenstates are expected to be of diffusive na- 
ture (in quasi-particle space), i.e. a natural generaliza- 
tion of the form postulated by Berry's conjecture. Such 
transition have been studied extensively in confined elec- 
tronic system, following a seminal paper by Althsuler, 



Kamenev, Levitov and Gefen (Altshuler et al. 1997) 



and on interacting electron systems with localized sin- 
gle particle states (Basko et al. 2006). While the de- 



pendence 77(A) was analyzed numerically in certain small 
systems (Canovi et al. 20ir;'Neuenhahn and Marquardt 



[2OIO; , Santos and Rigol^ 2010a,b ), eventual emergence of 



a sharp KAM-like threshold in the thermodynamic limit 
remains an open question. 

The ETH also suggested a new way on looking at quan- 
tum relaxational dynamics as dephasing in the many- 
body basis. In particular, the information about the 
asymptotic state is fully contained in the diagonal el- 
ements of the density matrix, which do not change in 
time if the Hamiltonian is constant. So the process of 



time evolution in this picture is equivalent to averaging 
of oscillating off-diagonal elements of the density matrix 
to zero. In a way this picture is different from conven- 
tional thinking based on kinetic theory of thermaliza- 
tion through collisions between quasi-particles and the 
time evolution of their distribution function. This ap- 
parent difference is hidden in the complicated structure 
of the many-body eigenstates. Our intuition is based on 
thinking about dynamics in the approximate basis, e.g. 
a basis of independent quasi-particles. The precise rela- 
tion between the many-body and kinetic approaches to 
thermalization is still an open question. Another poten- 
tially intriguing possibility is understanding thermaliza- 
tion as a renormalization group process, where time evo- 
lution results in averaging over high-energy degrees of 
freedom. If one deals with approximate noninteracting 
basis then because of interactions the process of elimi- 
nating high energy states affects dynamics of low energy 
modes and hence in renormalization of the low energy dy- 
namics. In Ref. (Mathey and Polkovnikov, 2010) it was 



shown that such renormalization process indeed can ex- 
plain real time dynamics in a two-dimensional sine Gor- 
don model and the emerging nonequilibrium Kosterlitz- 
Thouless transition. In Ref. ( Moeckel and Kehrein 2009 ) 



similar ideas were put forward to analyze dynamics of in- 
teracting fermions using the flow equation method. At 
the moment it is unclear whether using such real time 
renormalization group one can analyze relaxational long 
time dynamics in generic interacting systems. 



IV. EXPERIMENTAL PROGRESS IN QUANTUM 
DYNAMICS IN COLD ATOMS AND OTHER SYSTEMS 

As we mentioned before, the study of nonequilibrium 
dynamics of quantum many-body systems has been in- 
creasingly motivated by a series of advances in the field of 
ultracold atomic and molecular gases. Due to the conflu- 
ence of various features, these mesoscopic quantum sys- 
tems are in many ways near-ideal systems for the study 
of nonequilibrium quantum phenomena. 

Firstly, quantum gases can exhibit a remarkably high 
degree of isolation from environmental sources of deco- 
herence and dissipation. Thus, to an excellent approxi- 
mation, during duration of experiments they can be re- 
garded as closed quantum systems. Further, the dilute 
nature of these gases and exceptionally low temperatures 
result in long timescales of dynamical effects (typically 
on the order of milliseconds or longer) allowing for time- 
resolved studies of nonequilibrium processes resulting 
from phase-coherent many-body dynamics. Such stud- 
ies are hardly possible in conventional condensed matter 
systems. 

Secondly, an array of techniques have been developed 
to dynamically tune various parameters of the Hamil- 
tonian governing these quantum gases. This has made 
it possible to realize various prototypical nonequilibrium 
processes such as quantum quenches discussed above. 
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Quenches across phase transitions have been reahzed to 
investigate the onset and formation of long range order 
and the mechanism underlying the spontaneous forma- 
tion of topological defects. The latter is closely related to 
the KZ mechanism described earlier in the text. A quan- 
titative experimental study of the defect, entropy and en- 
ergy production resulting from such quantum quenches 
should allow for an accurate comparison with the theo- 
retical predictions. 

Lastly, the ability to engineer and experimentally real- 
ize model Hamiltonians of archetypal correlated systems 
coupled with a detailed knowledge of the microscopic in- 
teractions make ultracold atomic gases a tantalizing sys- 
tem for applications ranging from the quantum simula- 
tion of strongly correlated systems to the adiabatic quan- 
tum computation. In addition to the form of the model 
Hamiltonian, experimental control can also be achieved 
over the effective dimensionality of the ultracold gas mak- 
ing it possible to investigate the nontrivial interplay be- 
tween fluctuations, interactions and dimensionality. 

From a technological perspective, there is an increasing 
thrust towards engineering ultracold atomic many-body 
systems for applications in quantum metrology (Appel 



et al. 



et al. 



20091 lEstevelFal 



2008 Riedel et al. 



2008':'Lero ux et aH[20T0l|Meiser 
2010; Vengalattore ei aL||2007|. 



A deeper understanding of the dynamics of interacting 
many-body systems and the mechanisms of decoherence 
and dissipation in these systems is of crucial importance 
in this context. 

Motivated by these factors, a number of experiments 
have been performed in recent years using ultracold quan- 
tum gases to investigate topics including quantum coher- 
ent dynamics in optical lattices, quenches across quan- 
tum phase transitions and thermalization in low dimen- 
sional systems. For the purposes of this colloquium, 
we distinguish between classes of nonequilibrium exper- 
iments both in terms of the general protocol as well 
as the questions being addressed by these experiments. 
(i) Nonequilibrium states of many body atomic systems 
wherein the high degree of isolation of the atomic sys- 
tem from the environment allows for the creation of 
metastable or highly excited many-body states with long 
lifetimes, (ii) Quantum quench experiments in which 
one or more parameters of the Hamiltonian are changed 
rapidly to create an out-of-equilibrium state of the many- 
body system and (iii) Dynamical tuning of the Hamilto- 
nian in order to study quantum coherent dynamics of an 
interacting many-body system. 

These experimental advances have stimulated a very 
active theoretical research in the area of nonequilibrium 
quantum dynamics in interacting many-body systems. 
Among the issues most debated in recent literature is 
the relation between thermalization in isolated quantum 
systems and quantum integrability. In this regard, a re- 
cent pioneering study on thermalization in ID Bose gases 
was performed in Ref. (Kinoshita et al. 2006). In this 



atoms. The depth of the lattice potential far exceeded the 
energy of the ultracold gas ensuring negligible tunneling 
among the tubes. The array of tubes was then placed in a 
superposition of states of momentum ±2po by the appli- 
cation of a transient optical phase grating. The imparted 
kinetic energy was small compared to the energy required 
to excite the atoms to the higher transverse states and the 
gases remained one-dimensional. This out-of-equilibrium 
system was then allowed to evolve for variable durations 
before the momentum distribution was probed by absorp- 
tion imaging of the gas (see Fig. [2|. 

It was found that, while the initial momentum distri- 
butions exhibit some dephasing on account of trap an- 
harmonicities, the dephased distribution remains non- 
gaussian even after thousands of collisions. This is in 
distinct contrast to the gaussian distributions observed 
when the 2D optical lattice is adiabatically imposed on 
an equilibrium 3D Bose gas. This remarkable observation 
that the nonequilibrium Bose gases do not equilibrate on 
the timescales of the experiment appears consistent with 
the fact that this system is a very close experimental re- 
alization of a Lieb-Liniger gas with point-like coUisional 
interactions - an integrable quantum system in which 
only elastic pairwise interactions can occur. Apparently 
the experimental technicalities such as anharmonicities 
or the axial potential are insufficient to sufficiently lift 
integrability in this system. 

In addition to unambiguously showing the absence 
of thermalization within experimental timescales in this 
model realization of the Lieb-Liniger gas, this study also 
points the way towards addressing more general ques- 
tions on integrability and ergodicity. Starting from an 
integrable system, modifications such as the addition of 
finite range interactions, tunneling between the ID tubes 
or the imposition of axial potentials one can tunably lift 
integrability and analyze emergence of irreversability and 
thermalization. This experiment largely motivated much 
of the theoretical work discussed in the previous section. 

Another issue that has attracted a lot of attention is 
the search for universal effects either in the nonequilib- 
rium dynamics following a quantum quench or in the adi- 
abatic dynamics near a quantum critical point which we 
described earlier in the text. In particular, the issue of 
non-adiabatic dynamics near quantum phase transitions 
has been the focus of recent experimental studies on con- 
densate formation in a dilute, weakly interacting Bose 
gas that is rapidly c ooled past the BEC phase transition 
(Weiler et al. 2008). This process was found to be ac- 



experiment, a blue detuned 2D optical lattice was used to 
create arrays of tightly confined tubes of ultracold *^Rb 



companied by the spontaneous formation of topological 
defects, i.e. vortices, in the nascent superfluid. This can 
be phenomenologically understood as being due to the 
formation of isolated superfluid regions of a characteristic 
size 1^, each with a random relative phase. These isolated 
regions then gradually merge to give rise to global phase 
coherence. In this process, regions which enclose phase 
loops of 27r are constrained by the nature of the super- 
fluid, i.e. the continuity of the wavefunction, to have a 
vanishing superfluid density at the core. Thus, the KZ 
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FIG. 2 Time-of-flight absorption images of an ensem- 
ble of ID Bose gases - Ultracold atoms are confined in 
arrays of ID optical traps. Optical pulses are used to place 
the atoms in a superposition of ±2fifc momentum states. The 
gas is then allowed to evolve for variable durations before be- 
ing released from the trap and photographed to reveal the 
momentum distribution. The false color in each image is 
rescaled to show detail. The non-gaussian nature of the mo- 
mentum distribution clearly indicates an absence of thermal- 
iza,tion.( Adapted from (Kinoshita et al. 2006)) 



mechanism predicts a density of vortices that scales as 

i/e. 

In this experiment, a magnetically trapped thermal gas 
of ^''Rb atoms was cooled by radiofrequency (rf) evapo- 
ration to temperatures below the BEC transition tem- 
perature. The quench rate, i.e. the rate of cooling, was 
controlled by varying the rate at which the rf frequency 
was ramped down. Following a brief duration of equili- 
bration, vortices are detected by absorption imaging of 
the gas after ballistic expansion. Allowing for some un- 
certainty in the ability to discern a vortex due to line-of- 
sight integration in these images, it was found that about 
a quarter of the images showed at least one vortex core. 

The rate of cooling during the quench was limited by 
the collision rate between atoms in the trapped gas dur- 
ing evaporative cooling. This resulted in a limited dy- 
namic range for the quench rate. Also, the rapid de- 
crease of the thermal fraction following the formation of 
the condensate led to a low damping rate for the vortices. 
A faster quench rate, realized through a trap with tighter 
confinement or increased density or via sympathetic cool- 



ing with another species, could result in the observation 
of an increased number of vortices during the quench. In 
turn, this would potentially allow for quantitative tests of 
the predicted scaling of vortex number with the quench 
rate and the extraction of dynamic critical exponents. 

While the formation of a superfluid by quenching the 
temperature is seeded by thermal fluctuations, ultracold 
atomic gases also potentially allow for the realization of 
phase transitions initiated purely by quantum fluctua- 



tions ( Greiner et al. 2002a Sadler et al. 2006 1 



A particularly intriguing study of a quench past such 
quantum phase transitions was carri ed out in a degener - 



ate F = 1 spinor Bose gas of **^Rb (Sadler et al. 20061. 



These gases, with a spin degree of freedom arising from a 
non-zero hyperfine spin F , are quantum fluids that may 
simultaneously exhibit the phenomena of magnetism and 
superfluidity, both of which result from symmetry break- 
ing and long range order. Owing to rotational symmetry, 
the contact interactions between two atoms can be char- 
acterized by the total spin of the colliding pair. In the 
case of a i^ = 1 spinor gas, these interactions give rise 
to a mean field energy given by n(co + C2(F)^) where the 
coupling strengths co,2 are related to the s-wa ve scatter- 
ing lengths in the total spin / = 0, 2 channels (Ho 19981 



Ohmi and Machida 1998). In addition to the mean field 



interactions, a finite external magnetic field B imposes 
a quadratic Zeeman energy (QZE) that scales as q{F^) 
with q = (/XB_B)^/4A/,y where fiB is the Bohr magneton 
and Ahf is the energy splitting between the ground state 
hyperfine manifolds. 

For a F = 1 condensate of ^^Rb, the competing influ- 
ences of the spin-dependent interaction and the QZE give 
rise to a continuous quantum phase transition between a 
'polar' and a ferromagnetic phase. Rapidly tuning the 
external magnetic field from large values {q 3> |c2«|) to 
small values (q <^ \c2n\) quenches the spinor gas from 
the polar phase to the ferromagnetic phase. The ensu- 
ing growth of ferromagnetic domains was directly de- 
tected by in situ imaging (see Fig. [l]). It was found 
that the resulting texture of ferromagnetic domains was 
spatially inhomogeneous can characterized by a typical 
length scale that was related to the spin healing length 
^ — hj ^Jlm\c2n\. Concurrent with the appearance of 
these domains, the spin textures revealed the sponta- 
neous formation of polar-core spin vortices. These topo- 
logical defects are characterized by a non-zero spin cur- 
rent but no mass current. The origin of these spin vor- 
tices is also rooted in the KZ mechanism. It was shown 
that, for slow quenches, the number of such vortices is ex- 
pected to scale as r"^/^ where t is the time over which 



the spinor gas is swept into the ferromagnetic state ( Saito 



et al. 2007). 



The weak spin-dependent interactions inherent to this 
spinor gas also allow for nondestructive detection of the 
vortices and studies of their dynamics. In addition, the 
weak coupling between the spin and mass degrees of free- 
dom make it straightforward to realize extremely low 
spin temperatures to examine the role of quantum flue- 
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tuations in seeding this phase transition (Klempt et al. 
2010). These features make spinor quantum fluids a rich 



system to investigate the quench dynamics and KZ mech- 
anism past quantum phase transitions between different 
magneticaUy ordered phases. In addition, corrections to 



the KZ scaling imposed by long range interactions (Ven 



galattore et al. 2008 (, conservation laws and finite tem- 
perature effects can also be studied. 

Yet another range of experimental studies is made pos- 
sible by the tunability of atomic interactions using a Fes- 
hbach resonance. This technique allows the rapid dy- 
namic control of the s-wave scattering length by means 
of a time-varying external magnetic field. This ability 
was utilized in a recent study of a strongly interacting 
two-component Fermi mixture ( Jo et al. 2009 1 . Start- 



ing from an initially weak, repulsive interaction between 
the two Fermionic species, the interactions were rapidly 
increased by tuning the magnetic field to the vicinity of 
the Feshbach resonance. The subsequent decrease in the 
atomic loss rate, the increase in the size of the trapped 
gas and the increase in kinetic energy as measured in 
time-of-flight images were interpreted as an indication 
of the Stoner transition to a ferromagnetic state. How- 
ever, in a later theoretical work this interpretation was 
questioned and an alternative explanation based on rapid 



molec ule formation was suggested (Ref. (Babadi et al. 
2009)). Thus, a direct in-situ measurement of local mag- 



netization is necessary to understand whether ferromag- 
netism plays a role in this experiment. 

In addition to the thermalization dynamics across 
phase transitions, the long coherence times inherent to 
ultracold gases also makes it possible to study the quan- 
tum coherent dynamics of many-particle systems. A par- 
ticularly dramatic instance of such coherent many-body 
dynamics was illustrated in the collapse and revival of the 



matter wave field of a Bose condensate (Greiner et al. 



2002b). Here, the interaction-induced dynamical evolu- 



tion of a matter wave field was clearly revealed in the 
multiple matter wave interference patterns obtained af- 
ter releasing the gas from the lattice. This work has 
also been extended to the time-resolved observation of 
super exchange proce s ses in optical 'superlattice' poten- 
tials (Trotzky et al. 20081. Similar demonstrations of 



coUisional coherence have also been shown in spinor Bose 
gases (Chang et al. 2005 Kronjager et al. 2005). Due 



to the internal degrees of freedom in a spinor gas, the 
dynamics in this fluid is due to coherent spin-mixing col- 
lisions. In a trapped gas that is well described by the 
single mode approximation (SMA), these coherent colli- 
sions can lead to the periodic and reversible formation of 
condensates in initially unpopulated spin states. 

Further, in certain situations, coherent interactions can 
also lead to quantum correlations (Sorensen et al. 2001). 



Schemes that might realize such entangled many-particle 
states have received attention due to potential applica- 
tions in quantum information processing and metrology. 
The dynamical evolution of such entangled states in the 
presence of quantum or thermal fluctuations is obviously 



of great interest. A recent experiment investigated this 
evolution in low dimensional two-component Bose gases 



with adjustable interactions ( Widera et al. 



20081, finding 



that quantum fluctuations play a crucial role in the phase 
diffusion dynamics of low dimensional systems. More 
recently, the dynamical control of a Bose-Einstein con- 
densate confined in a strongly driven optical lattice was 
demonstrated ( jLignier et al.\ 2007 ) . By periodically mod- 
ulating the lattice potential, the tunneling parameter J 
was shown to be suppressed in a phase coherent manner 
opening the possibility of driving quantum phase transi- 
tions using this technique. 

The isolation of ultracold atomic gases from external 
sources of dissipation also makes it possible to study re- 
laxation dynamics driven purely by intrinsic mechanisms. 
Such mechanisms should set the timescales for adiabatic 
quantum computing or the simulation of strongly cor- 
related lattice models. A recent experiment along these 
lines investigated the evolution of excited states of the re- 



pulsive Fermi-Hubbard system (Strohmaier et al. 20101. 



Here, doubly occupied lattice sites (doublons) were cre- 
ated by modulating the lattice and the subsequent decay 
of the system to thermal equilibrium was monitored over 
time. It was shown that the lifetime of these doublons 
scales exponentially with the ratio of the interaction en- 
ergy to kinetic energy, in fair agreement with theoretical 
predictions. It was argued that the dominant mechanism 
driving this relaxation was a high-order scattering pro- 



cess involving several fermions (Sensarma et al. 2010) 



While this colloquium places an emphasis on experi- 
ments involving ultracold atomic gases, there is a range of 
other mesoscopic quantum systems which also lend them- 
selves to studies on quantum nonequilibrium dynamics. 
For the sake of completeness, we briefly review a few of 
these systems here. Defect formation following a quench 
was first studied in the context of vortices in liquid crys- 
tals (Chuang et al. 1991). This has since been followed 



by similar studies in various mesoscopic systems includ- 
ing isolated superconducting loops where the defects as- 
sume the form of spontaneous fiuxoids ([Monaco et al. 



and multi-Joscphson junction loops (Monaco et al. 12006 



2009), superconducting thin films (Maniv et al.[ |2003| 



2002). A cumulative view of these studies would indicate 
that the infiuence of finite size effects, thermal fiuctua- 
tions and dimensionality on the production of topological 
defects by the KZ mechanism is as yet unclear and a topic 
that warrants further study. 

Another potential system for the study of nonequilib- 
rium dynamics of many-particle states arises from rapid 
advances in the field of photonics. There have been sev- 



eral proposals ( Chang et al. 2008 Greentree et al. 2006 ) 



for the dynamical creation of strongly correlated pho- 
tonic states using photon-photon interactions mediated 
by a nonlinear optical medium. The realization of states 
such as a Tonks gas of photons have been proposed using 
hollow-core optical fibers, tapered optical fibers, photons 
in coupled cavities and surface plasmons on conducting 
nanowires. Such strongly correlated photon states should 
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have applications in metrology, sub-shot noise interfer- 
ometry and the quantum emulation of exotic spin mod- 
els. 



V. OUTLOOK 

One of the ultimate goals of the new field of quantum 
dynamics is to develop a systematic understanding of 
nonequilibrium phenomena in strongly interacting quan- 
tum many-body systems. A few of the most significant 
open questions along this avenue are readily identified : 
How can we classify nonequilibrium behavior in closed 
many-body systems? What is the general relation be- 
tween integrability and dynamics? What is the dynam- 
ical effect of a gradual breaking of integrability? What 
are the effects of dissipation on these nonequilibrium pro- 
cesses? Can we understand time evolution of interacting 
systems through the renormalization group? Answering 
these and other questions allied with systematic, quanti- 
tative studies of possible nonequilibrium quantum phase 
transitions and the extraction of dynamical critical ex- 
ponents, are just a few of the many tantalizing programs 
to be pursued. The rapidly developing sophistication 
and precision of ultracold atomic experiments and other 
experimental systems should allow for close and direct 
comparison between theoretical predictions and ad hoc 
experiments. 

The realization of robust techniques for the experimen- 
tal study of such systems and the development of theo- 
retical tools to describe nonequilibrium many-body pro- 
cesses should bode for tantalizing opportunities in this 
nascent field, potentially leading to a deeper understand- 
ing of the principles governing nonequilibrium many- 
body phenomena and establishing robust connections be- 
tween microscopic dynamics and statistical physics. 
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